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As  part  of  the  on-going  program  in  "Decision  and  Control  Models  in 
Operations  Research,"  Mr.  Ronald  Kaufman  investigates  the  behavior  of 
multi-item  inventory  control  systems  in  which  the  means  of  the  under- 
lying demand  distributions  vary  in  a cyclic  manner,  corresponding  to 
seasonal  peaks  for  inventoried  goods.  Mr.  Kaufman  adapts  the  Ehrhardt 
power  approximation  (Technical  Report  No.  7),  which  was  originally  de- 
signed for  stationary  demand  distributions,  to  the  nonstationary 
environment.  Mr.  Kaufman  tests  the  behavior  of  the  power  approximation 
in  inventory  systems  where  the  demand  distribution  parameters  are 
estimated  from  a limited  amount  of  historical  information.  Several 
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Other  related  reports  dealing  with  the  research  program  are  given  below. 
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ABSTRACT 


(s,S)  INVENTORY  POLICIES  IN  A 
NONSTATIONARY  DEMAND  ENVIRONMENT 

Ronald  Louis  Kaufman 
University  of  North  Carolina-1977 

Scarf  (1960)  proves  the  optimality  of  (s,S)  policies  for  a class 
of  discrete  review  nonstationary  inventory  models.  A considerable 
amount  of  inventory  literature  concerns  computation  of  optimal  and  ap- 
proximately optimal  (s,S)  policies  under  the  Scarf  hypotheses.  Little 
research  has  dealt,  however,  with  the  case  of  nonstationary  demands. 

This  investigation  examines  the  situation  in  which  demand  distributions 
are  independent,  but  not  identically  distributed,  and  vary  in  a cyclic 
manner.  Products  that  experience  seasonal  demands  are  a typical  example 
of  such  a demand  process. 

A detailed  analysis  of  the  nature  of  optimal  policies  in  the  non- 
statlonary  environment  is  presented.  The  behavior  of  selected  operating 
characteristics  such  as  period-end  inventory,  backlogged  demand,  fre- 
quency of  stockout,  replenishment  frequency,  and  associated  costs  also 
are  examined.  The  above  performance  measures  for  an  individual  inven- 
toried good  are  aggregated  over  all  goods  to  provide  an  analysis  of  the 
multi-item  system  behavior. 

Approximately  optimal  (s,S)  policies  are  derived  for  the  nonstation- 
ary environment.  The  policies  are  based  on  the  power  approximation  of 
Ehrhardt  (1976),  and  require  knowledge  of  only  the  mean  and  variance  of 
demand.  The  operating  characteristics  of  approximately  optimal  policies 
are  compared  with  those  of  optimal  policies. 


The  approximately  optimal  policy  rule  is  examined  in  a statistical 
environment  that  generalizes  MacCormick  (1974).  Policy  parameters  are 


revised  periodically  using  a limited  history  of  past  demands  to  estimate 
the  mean  and  variance  of  demand.  Each  time  the  policy  parameters  are 
revised,  forecasts  of  system  operating  characteristics  are  calculated 
from  a retrospective  simulation  employing  the  same  sample  of  demands 
that  was  used  to  revise  the  policy. 

The  statistical  phenomena  are  studied  by  means  of  a computer 
simulation  program  using  time  series  analysis  techniques. 
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1.  INVENTORY  MANAGEMENT  WITH  INCOMPLETE 
DEMAND  INFORMATION 


This  section  describes  the  inventory  model  assumptions  and  the 
inventory  operating  characteristics  that  will  be  used  to  evaluate  in- 
ventory policies.  The  Ehrhax  t power  approximation  is  also  presented, 
and  we  will  show  how  it  can  be  successfully  adapted  to  inventory  systems 
with  nonstationary  demands. 

1.1  The  Model 

1.1.1  Structure 

We  consider  a single-item  inventory  model  in  which  the  inventory 
level  is  reviewed  each  period,  and  demand  for  the  item  is  generated 
by  a discrete  time  stochastic  process.  The  demand  process  is  non- 
stationary; Specifically,  the  demands  are  independent,  but  not  identi- 
cally distributed,  non-negative  integer-valued  random  variables.  The 
distributions  used  to  generate  these  demands  are  cyclic  with  period  L. 
Let 


51*  *>2"“’  ^L*  ^L+l  * * * * * ^2L*  f’2L+l’ 


•,C3L*  ^ 3L+1 * ' ' ‘ 


represent  the  demand  sequence:  then  the  joint  distributions  of  CnL+1t 

e(n+l)L  and  Wl’-"’  ^(m+l)L  are  are  a8Sumed  t0  be  the  same  for 
any  integers  m and  n _>  0.  Since  there  are  at  most  L different 

demand  distributions,  4>^,  for  i » 1,...,L,  denotes  the  distribution 

function  for  the  i*’*1  demand  within  the  cycle.  We  assume  that  demand 


is  met  as  long  as  there  is  stock  on  hand,  and  when  a stockout  occurs, 
unfilled  demand  is  backlogged  until  sufficient  replenishment  arrives. 
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Units  of  inventory  are  conserved,  there  being  no  losses  by  de- 
terioration, obsolescence,  or  pilferage,  and  disposal  is  not  allowed. 
Inventory  on  hand  at  the  end  of  a current  period  is  the  inventory  from 
the  previous  period  plus  any  replenishment  that  a-rives,  less  demand 
in  the  current  period.  If  inventory  on  hand  is  negative,  the  amount 
represents  backlogged  demand.  Replenishments  are  assumed  to  be  de- 
livered a fixed  leadtime  X after  being  ordered.  The  time  sequence  of 
events  within  any  period  is  taken  to  be  order,  delivery,  and  de- 
mand. 

The  cost  functions  for  the  model  are  simple.  With  no  discount- 
ing over  time  and  an  unbounded  planning  horizon,  a resonable  criter- 
ion is  minimization  of  long-run  average  expected  total  cost  per  period. 
The  cost  of  a replenishment  quantity  is  assumed  to  be  linear  with  a 
fixed  ordering  cost  K.  The  inventory  holding  cost  is  proportional  to 
any  stock  on  hand,  at  unit  cost  h,  assessed  at  the  end  of  each  period; 
a unit  penalty  cost  it  is  applied  to  any  quantity  on  backorder  at  the 
end  of  each  period. 

The  associated  control  problem  is  to  make  optimal  replenishment 
decisions  for  this  model.  We  assume  that  control  over  replenishment 
is  exercised  by  an  (s^,S^)  policy,  for  i - 1,...,L:  Whenever  in- 
ventory x on  hand  and  on  order  at  the  start  of  the  i^  period,  1 
i £ L,  within  a cycle  of  length  L,  drops  below  the  value  s^,  we 
place  a replenishment  order  of  size  S^-x.  When  the  objective  of  con- 
trol is  to  minimize  Jong-iun  average  expected  total  coat  per  period, 
and  the  distributions  and  parameters  specifying  the  structure  of  the 
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model  are  completely  known,  it  follows  from  the  theory  of  dynamic  pro- 
gramming, [see  Veinott  (1966)  and  Karlin  (1959)],  that  there  is  an 
(Sj,S^)  policy,  for  i - 1,...,L,  that  is  optimal  among  all  possible 
policies . 

When  the  distribution  generating  demand  is  not  know,  then 
even  though  this  is  the  only  assumption  relaxed,  it  may  no  longer 
be  true  that  a (s^,S^),  for  i * 1,...,L,  policy  is  optimal,  using 
the  criteria  of  long-run  average  expected  total  cost  per  period. 
Nevertheless,  in  this  study  we  shall  use  3uch  a policy,  since  it 
is  in  popular  use  in  the  applied  situation  of  incomplete  informa- 
tion. Information  about  the  L demand  distributions  4^,  for 
i = 1,...,L,  must  be  gathered  by  inference  from  a limited  sample 

of  demands. 

1.1.2  Algorithms  to  Set  (s,S)  with  Partial  Knowledge 

2 

When  the  means  and  variances  o^,  for  i = 1,...,L,  of 

each  demand  distribution  in  a cycle  are  known,  approximations  to 
the  optimal  policy  are  available.  The  approximations  are  based 
on  asymptotic  theory  that  assume  large  unit  penalty  costs  and  large 
replenishment  costs.  The  algorithm  adopted  here  is  the  power 
approximation  of  Ehrhardt  (1976),  which  is  a generalization  of 
earlier  policies  based  on  the  work  cf  Roberts  (1962).  Ehrhardt  extend- 
ed Roberts'  work  by  generalizing  the  functional  forms  for  s and  S-s, 


4 


4 


and  then  fitting  these  functional  forms  to  a larger  set  of  known  optimal 
policies  using  least-squares  regression.  The  power  approximation  was 
designed  by  Ehrhardt  for  use  in  stationary  environments;  one  of  the 
purposes  of  this  study  is  to  investigate  adaptations  of  the  power 
approximation  to  nonstationary  environments. 

Let  p^X  + ^and  for  i “ 1,...,L,  denote  the  arithmetic 

average  of  the  next  X+l  demand  means  and  variances,  respectively, 
starting  in  the  i1*1  period  of  the  cycle.  For  example,  if  i - 5 and 


X = 2,  then  p 


(3) 

5 


-2(3) 

0 


and  oP"  are  the  arithmetic  averages  of  the  means 


and  variances,  respectively,  of  the  demands  in  periods  5,  6,  and  7.  For 

—2 (X+l) 


, . , — (X+l) 

each  pair  of  p^,  and 


the  power  approximation  computes. 


(1) 

and 

(2) 

where 


d . «.*63)(r;i+i))-364  (K/h)-498to+i)«3(1+1)]-°691 


i 

(X+l)  p<X+1)  + [(X+Dp^']™  (a;'Ari//ur,i,)-''''vU(z) 


,-(X+l)-,.416  ,-2(X+l)-  (X+l),. 603, 
i v i 


Si  ” 8i  + Dp  * 


U(z)  is  given  by 


U (z)  - (. 182/z)  + 1.142  - 3.466z, 


(3) 


(^- (X+l))  .364  (K/h).498 
(1  + [(X+l)  o2<A+1)]-431 


r (X+l) 


If  D i s greater  than  1.5,  let  s^  and 


Otherwise,  compute 
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(M-l)u 


(X+l) 


-2(X+l)-,»s 


+ v [(x+1) 


where  v is  the  solution  to 


F (v)  * u/(u+h)  , 


and  F(*)  is  the  cumulative  distribution  function  of  the  unit  normal 
distribution.  The  policy  parameters  are  then  given  by 
s.  = minimum  {s^  , S^} 

(6) 

S = minimum  {S^  , S^} 

Since  demands  are  integer-valued,  we  round  s^,  D^,  and  to  the 

nearest  integer. 

The  policy  (s^.S^)  is  utilized  at  the  beginning  of  the  i1^  period  in 
the  cycle.  The  policy  (s^.S^),  for  i ■ 1,...,L,  computed  using  the  power 
approximation,  has  been  adopted  in  this  study  as  the  basic  decision  rule 
in  the  case  of  limited  information  about  the  demand  distribution  para- 
meters . 

1.2  Experimental  Design 
1.2.1  Nonstationary  Demands 

In  this  study,  we  assume  that  an  item's  demand  distribution  in  any 
period  is  negative  binomial.  This  distribution  has  parameters  r and 
p,  where  r>0  and  0<p<l,  with  probability  mass  function 


f(x)  * ro£)Xr(r)  U-p*"  pX  for  x “ °>l 

yielding  mean  and  variance  — ^ — 5-  . 

1_P  (1-p) 

Note  that  the  var iance-to-mean  ratio  is  -jpy  , which  is  a function  of  p , 
and  throughout  this  study  the  variance-to-mean  ratio  is  fixed  at  3;  thus 
p ■ 2/3.  We  introduce  nonstationarity  by  varying  the  parameter  r 
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cyclically.  We  select  cycle  length  L *»  12,  corresponding  to  a yearly 
demand  cycle  of  12  months'  demands. 

We  select  for  study  five  general  types  of  cyclic  demand  structures. 
Figure  1.1  illustrates  the  movement  of  the  demand  distributions'  means 
in  each  cycle.  Figure  1.1b  displays  a step  model  demand  structure  that 
henceforth  we  refer  to  as  model  I.  The  mean  demand  level  for  this  model 
is  a for  the  first  nine  periods  of  the  cycle,  and  2a  for  the  last 
three  periods  of  the  cycle. 

A more  gradual  increasing  mean  is  shown  for  models  II  and  III  in 
Figures  1.1. c and  1.1. d,  respectively.  Both  models  maintain  a constant 
mean  a for  the  first  seven  periods  of  the  cycle.  The  mean  gradually 
increases  to  a maximum  level  in  period  10  of  2a  for  model  II  and  5a 
for  model  III.  After  period  10  the  mean  demand  slowly  decreases. 

Model  IV,  depicted  in  Figure  1.1. e shows  mean  demand  varying  in  a 
sinusoidal  manner  where  a is  the  lowest  level  of  mean  demand  and  5a 
is  the  largest  level. 

Throughout  this  study  models  I,  II,  III,  and  IV  are  compared  with 
the  stationary  model,  depicted  in  Figure  1.1. a.  The  stationary  model 
maintains  a constant  mean  demand  level  a . The  five  models'  values 
of  are  selected  to  ensure  that  overall  mean  period  demand  for  the 
entire  cycle  is  equal  to  values  specified  in  Section  1.2.2. 

1.2.2  Design  Parameters 

The  single-item  model  described  in  Section  1.1.1  is  studied  under 
the  range  of  input  parameters  given  in  Table  1.1.  The  demand  distribu- 
tion is  negative  binomial  with  variance-to-mean  ratio  fixed  at  3.  The 
two  mean  values  for  demand  are  8 and  16.  Two  values,  2 and  4,  are 
assigned  to  leadtime.  Since  the  cost  function  is  linear  in  the  para- 
meters K,h,  and  *,  the  value  of  the  unit  holding  cost  is  set  at  unity. 
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Figure  1.1  Demand  Means 


Essential  changes  in  cost  arise  only  with  changes  in  the  ratios  K/h 
and  ir/h.  The  stockout  costs  are  n = 4 and  99,  and  the  fixed  cost 
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values  are  K = 32  and  64.  All  combinations  of  input  parameters  are  ex- 
amined in  the  study. 

Each  set  of  16  input  combinations  is  examined  under  each  of  the 
five  demand  structures  discussed  in  Section  1.2.1,  resulting  in  five 
separate  inventory  systems.  The  systems  are  labeled  model  I,  model  II, 
model  III,  model  IV,  and  the  stationary  model,  corresponding  to  the 
names  given  to  the  demand  structures.  These  five  systems  constitute 
our  basic  experiment. 


Table  1.1 
System  Parameters 


Parameter 

Parameter  Settings 

Number 
of  Settings 

Demand  distribution 

Negative  Binomial  (V/y  - 3) 

1 

Mean  demand 

8,  16 

2 

Unit  holding  cost 

1 

1 

Unit  backlog  penalty  cost 

4,  99 

2 

Replenishment  setup  cost 

32,  64 

2 

Replenishment  leadtime 

2,  4 

2 

1.2.3  Decision  Rules 

This  study  examines  three  decision  rules:  optimal  inventory  rules, 
approximately  optimal  rules,  and  statistical  approximation  rules.  The 
computation  of  optimal  policies  is  explained  in  detail  in  Section  2.1. 
Approximately  optimal  policy  rules  are  calculated  using  the  power 
approximation,  discussed  previously  in  Section  1.1.2.  A third  method 
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is  the  statistical  power  approximation,  which  is  explained  in  the 
remainder  of  this  section.  Both  the  optimal  and  approximately  optimal 
policies  are  analyzed  analytically  in  Section  2.  Section  2 also  inves- 
tigates the  effects  of  demand  structure  misspecif ication  on  the  power 
approximation.  Computer  simulation  is  used  to  investigate  the  statistical 
power  approximation,  and  the  findings  are  presented  in  Section  3. 

For  an  arbitrary  integer  t,  let  G<t  be  the  number  of  observed 
demands  during  periods  t-G  to  t-1.  In  order  to  use  the  statistical  power 
approximation  of  Ehrhardt  (1976)  in  a stationary  environment,  the  sta- 
tistics required  are  the  sample  mean  and  variance  of  demand,  C and  V, 
respectively.  If  t is  the  period  in  which  we  decide  to  compute  a 
statistical  policy,  then 

- -1  G 

C “ G EC 

T-1  t-T 

-1  G - 2 
V - (G-l)  1 E (C-  O . 

T-1 

When  using  the  statistical  power  approximation,  we  periodically  obtain 
a new  (s,S)  policy  by  substituting  C and  V for  and 

in  equations  (1)  through  (6). 

In  a regression  context  C and  V are  estimators  of  the  regression 
coefficient  a and  the  variance  of  the  disturbance  term  for  the  re- 
gression model 


y - ox  + e for  t - 1,...,G, 
•'t-T  t-T  t-T  ’ ’ * 


where,  for  the  stationary  environment 


yfc  T • demand  in  period  t-T 


L 
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a m stationary  mean  demand  level 


x - 1 
t-T 


E(et_T)  = 0 


for  all  t>T 


E(e  e^_  ) * 0 for  t^k,  t>T,  k>T, 


E(Vt)  “ °2 


for  all  t>T  . 


In  our  nonstationary  environment,  however,  x£_T  Is  no  longer 
Identically  1 for  every  demand  period;  now  xt_T  depends  on  the  posi- 
tion of  period  t-T  within  the  demand  cycle.  We  assume  that  we  know 
the  cycle  length,  and  we  let  each  period's  mean  demand  In  the  cycle  be 
represented  as  the  product  of  the  appropriate  value  of  x^  ^ and  a , 
which  is  the  smallest  demand  mean  of  the  demand  structure.  We  also 
assume  that  the  values  of  x^  are  known,  but  that  a must  be  estimated. 
The  appropriate  regression  model  is 


y = otx^  + e 
t-T  t-T  t-T 


for  t - 1, . . . T ( < t ) , 


where 

y T - demand  In  period  t-T 

a =■  smallest  mean  demand  level  In  the  demand  structure 

x.  - appropriate  constant  such  that  ooc  is  the  mean 
t-T  r t-T 

demand  level  in  period  t-T 

E(et  t ) • 0 for  all  t>T, 

E(et_Tek_T)  - 0 for  t+k,  t>T,  k>T, 

2 2 

E(c ) - ^(t-t)  for  all  t>T, 


- c for  all  t>T , 


y(t-T)  - appropriate  demand  structure  index  used  in  period  t-T 
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y (t-t)  Identifies  the  demand  distribution  used  in  period  t-x,  relative  to 
the  cycle.  For  example,  if  y(t-x)  - 5,  then  period  t-x  receives  demand 

from  the  fifth  period  of  the  cycle,  and  the  variance  of  demand  in  period 
. 2 

t-T  iS  Oy 


Table  1.2 
Values  for  ^ 

Period  In  Cycle 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Stationary  Model 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Model  I 

1 

1 

1 

1 

1 

1 

1 

1 

X 

2 

2 

2 

Model  II 

1 

1 

1 

1 

1 

1 

1 

1.33 

1.67 

2 

1.67 

1.33 

Model  III 

1 

1 

1 

1 

1 

1 

1 

2.33 

3.67 

5 

3.67 

2.33 

Model  IV 

2.60 

1.80 

1 

1 

1.80 

2,60 

3.40 

4.20 

5 

3 

4.20 

3.40 

The  values  of  xt  within  a cycle  for  each  of  the  demand  structures 

studied  are  given  in  Table  1.2.  Observe  that  x^  is  1 when  the  mean  demand 
level  in  period  t-x  is  a and  2 or  5 for  the  nonstationary  systems  when 
the  mean  demand  level  in  period  t-x  is  at  its  largest  value.  Hence,  in 
models  I through  IV,  a represents  a base  level  of  mean  demand,  and 
determines  how  much  above  that  level  the  mean  demand  in  period  t-x  will  be. 

i 

The  nonstationary  regression  model  contains  heteroscadastic  disturb- 
ances. However,  dividing  both  sides  of  the  regression  equation  by  /x^  ~ 
will  eliminate  the  heteroscedasticity  [see  Johnston  (1972)],  and  produce 
a transformed  model  with  disturbance  variance  equal  to  the  variance-to- 
mean  ratio,  which  has  been  kept  constant,  times  a . The  resulting 
least  squares  estimator  of  a is 


G 

E y 

T-l 


t-x 


G 
E K 
x-1 


t-x 
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and  a standard  estimator  of  ca  the  variance- to-mean  ratio  times  a , is 

£ ~ (yc-r  - xt-T)2 

2 T-l  t-T  C T tT 


G - 1 

Once  a and  s are  obtained,  we  can  estimate  and  cr^  , 

for  i - 1, . . . ,L,  by 

R+X^ 

£ a x, 

5^+1)  , L 

wi  X + 1 


and 


, s2-(»+l) 


1 ' ~l 

respectively,  where  R is  chosen  such  that  Y(R)  * i.  These  estimates  are 

substituted  for  and  a^+1^  in  equations  (1)  through  (6)  of 

Section  1.1.2  to  obtain  the  statistical  policy  (s^,S^). 


1.2.4  Policy  Revision 

In  implementing  a statistical  inventory  policy,  in  a real-life 
environment,  the  policy  parameters  s^  and  for  i ■ 1,...,L,  would  be 
revised  at  regular  intervals.  With  such  revisions  taking  place,  there 
are  system  design  choices  to  be  made  concerning  the  length  of  history 
kept  and  the  amount  of  information  stored.  For  the  simulation  experi- 
ments in  this  study,  a history  of  fixed  length  is  kept  for  each  re- 
vision, and  equal  weight  is  given  to  each  observation.  This  is  not 
optimal  if  the  demand  process  is  known  to  be  stationary,  for  then  the 
entire  history  should  be  accumulated  to  give  progressively  better  know- 
ledge and  performance.  Even  when  demand  is  known  to  be  nons Cat ionary,  but 
varying  in  a regular  manner,  such  as  by  a trend  or  periodic  cycle,  or  both, 
an  optimal  decision  rule  would  generally  utilise  the  entire  history. 

u 

•• I ■ I— ■ - — ■ — — ■ ■ ■ — — 
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The  decision-maker  usually  Is  not  in  a position  to  know,  however, 
that  conditions  observed,  even  over  the  entire  past  history,  will  con- 
tinue to  prevail.  This  lack  of  complete  information  provides  justifi- 
cation for  making  frequent  revisions,  placing  greater  weight  on  observa- 
tions from  the  immediate  past  and  less  on  distant  history.  For  this 
study,  an  arbitrary  choice  has  been  made  to  keep  a history  of  fixed 
length  and  give  equal  weight  to  all  observations  in  the  history.  In 
this  study,  the  revision  interval  and  revision  history  length  are  both 
24  periods,  implying  "monthly"  demand  data  that  are  accumulated  over  the 
same  time  interval.  For  a detailed  discussion  of  demand  generation, 
system  initialization,  and  output  analysis,  in  the  simulation  experiment 
see  MacCormick  (1974,  pp.  52-58). 

1.2.5  Item  Operating  Characteristics 

We  examine  detailed  operating  characteristics  of  each  item  in  the 
five  systems.  These  are  the  period-end  inventory  on  hand,  the  period- 
end  stockout  quantity,  the  frequency  of  period-end  stockouts,  the  re- 
plenishment quantity,  the  frequency  of  replenishment,  and  the  total  cost 
incurred.  It  is  assumed  that  the  decision-maker  is  interested  only  in 
the  average  performance  of  these  characteristics  between  revisions, 
corresponding  to  the  practice  adopted  by  accountants  of  making  reports 
to  management  at  periodic  intervals. 

One  objective  of  the  study  is  to  make  inferences  about  the  dis- 
tributions of  these  operating  characteristics  for  each  of  the  three 
decisions  rules  (optimal  policies  under  full  information,  approximately 
optimal  policies,  and  statistical  approximation  policies),  that  were 


discussed  in  Section  1.2.3. 
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1.2.6  Forecasting 

The  method  of  MacCormick  (1974)  is  used  to  forecast  operating  char- 
acteristics when  the  system  is  controlled  by  the  statistical  power  approx- 
imation rule.  Often  termed  "retrospective  simulation,"  this  method  takes 
the  recent  demand  history  for  the  items  in  the  sample  and  estimates  the 
performance  of  the  chosen  policy  for  this  history.  Thus,  the  method  uses 
the  same  data  twice,  once  to  fix  the  policy  parameters,  and  once  to  fore- 
cast the  performance.  As  a result,  the  forecasts  are  biased,  and  this 
study  investigates  the  extent  of  the  bias. 

At  each  revision,  after  the  (s^,S^)  parameters,  for  i = i,...,L, 
have  been  set  to  new  values,  a forecast  is  made  of  the  properties  of  the 
system.  A history  of  the  latest  g+*  demand  values  is  kept  to  make  the 
forecast.  [We  note  that  there  are  differences  in  information  storage 
requirements  of  considerable  importance  in  practice.  To  make  a forecast 
of  performance  by  retrospective  simulation,  the  demands  for  the  item  In 
each  of  G+*  periods  must  be  stored,  whereas  the  statistical  decision  rule 
may  require  only  a handful  of  sufficient  statistics  to  set  the  (s^,S.) 
policies.  Also,  for  multi-item  systems  the  cost  of  accumulating  in  stor- 
age the  data  needed  to  make  forecasts  by  retrospective  simulation  is 
likely  to  be  so  high  that  histories  will  be  kept  for  only  a representative 
sample  of  items  from  the  system.] 

The  forecasts  are  made  every  24  periods  by  running  the  system,  using 
the  new  policy  values  of  s^  and  , for  i ■ 1,...,L,  and  observing 
the  operating  characteristics  as  if  these  policy  values  had  been  in  force 
when  the  history  occurred.  Each  time  there  is  a policy  revision,  (every 
24  periods),  there  is  a forecast  for  each  operating  characteristic.  The 
simulation  initializes  the  stock  on  hand  and  on  order  at  the  actual  value 


at  the  time  of  the  forecast,  that  is,  at  the  time  of  revision.  As  actual 
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stock  on  hand  is  not  recorded  until  the  elapse  of  a leadtime  after  re- 
vision, the  same  interval  is  allowed  to  elapse  before  recording  it  for 
the  forecast.  The  inventory  on  hand  variable  is  therefore  initialized 
at  the  initial  value  for  inventory  on  hand  and  on  order  less  the  first 
X demands  In  the  history  of  G+X  demands. 

1 . 3 A Multi- item  System 

Scientific  techniques  for  inventory  control  are  generally  applied  to 
systems  of  many  items.  This  study  combines  the  results  from  the  five  16 
single-item  systems  into  five  multi-item  systems.  Since  management  general- 
ly assesses  the  performance  of  control  techniques  for  a multi-item  inven- 
tory system  by  observing  indices  that  are  aggregate  operating  characteris- 
tics [Wagner  (1962)3  » certain  aggregate  characteristics  have  been  computed. 

The  operating  characteristics  of  the  multi-item  system  under  statis- 
tical control  have  been  measured  by  aggregating  the  sample  values  of  the 
corresponding  characteristics  for  each  item  in  the  system.  When  the 
system  is  operated  under  perfect  information,  these  characteristics  are 
computed  analytically.  The  aggregate  of  average  total  cost  per  period  Is 
computed  as  the  arithmetic  sum  of  the  corresponding  costs  for  each  item. 

The  components  of  total  cost  for  inventory  storage,  backlog  penalty,  and 
ordering  replenishments  are  similarly  computed.  The  aggregate  backlog  and 
replenishment  frequencies  are  arithmetic  averages  of  the  corresponding 
frequencies  observed  for  each  item  in  the  system.  Since  the  unit  Inven- 
tory holding  cost  for  all  items  is  unity,  the  average  number  of  units  In 
inventory  at  period-end  is  numerically  identical  to  the  aggregate  average 
holding  cost  per  period.  Finally,  a weighted  proportion  of  demand  back- 
logged  is  computed  as  the  ratio  of  a weighted  sum  of  the  average  quantity 
backlogged  per  period  to  a weighted  sum  of  the  (exact)  mean  values  of 
demand.  The  weights  used  in  both  the  numerator  and  denominator  of  the 
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ratio  are  the  unit  cost  of  backlogging  demand  for  the  respective  Item. 
Throughout  this  study  the  nonstationary  systems  will  be  compared  to 
the  system  with  stationary  demands. 
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2.  SYSTEM  CONTROL  WITH  FULL  INFORMATION  AUOUT  DEMAND 

Section  2.1  contains  a detailed  explanation  of  the  analytical  models 
used  to  obtain  optimal  policies  and  evaluate  the  operating  characteristics 
of  arbitrary  policies  of  the  (s,S)  type.  Sections  2.2  and  2.3  record 
information  descriptive  of  single-item  and  multi-item  inventory  systems 
controlled  optimally  and  approximately  optimally  with  full  information 
about  demand.  The  main  intention  of  Section  2 is  to  establish  bench- 
mark values  of  the  systems'  operating  characteristics  foi  later  comparison 
with  systems  controlled  using  only  statistical  Information  about  demand. 
The  data  in  this  section  also  illustrates  that  the  Ehrhardt  power  approx- 
imation can  be  adapted  to  control  nonstationary  inventory  systems. 

2 . 1 Analytical  Models  Used  to  Obtain  and  Evaluate  Optimal  Policies 

In  Section  1.1.1  we  denoted  the  length  of  the  demand  distribution 
cycle  by  L.  Hence  If  represents  a demand  density  for  the  i 

period  in  a cycle  of  L periods,  then  demand  densities  in  successive 
periods  are 

♦ $2’  • * * * ♦p  4>2»  • • • . ♦j  » • • • • 


Thus,  is  the  demand  density  in  period  1,  $ is  the  demand 

density  in  period  2,  and  is  the  demand  density  in  period  L. 

In  period  L+l,  the  demand  density  is  again  <f> ^ , and  the  cycle  repeats. 

We  denote  the  v+1-  fold  convolution  of  successive  demand  densities  by 
v+1  v+1 

<f> ^ ; that  is,  is  the  convolution  of  i and  the  next  v densities 

that  follow  i in  the  cycle,  for  i - 1,...,L. 

When  n periods  remain  in  the  finite  horizon  of  a single-item 
inventory  model,  an  optimal  (s^.S^)  policy  is  obtained  by  employing  a 
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recursive  relationship  based  on  the  results  of  Scarf  (1960)  and  Veinott 
(1966).  Wt*  utilize  the  relationship  for  n > A 


f (x)  = min  {K6(y  - x)  + H(y,y(n))  + 5.  f (y  - £)<t>  , A O } 

n y^x  C=0  n~l  Y(n) 

where 

x = inventory  on  hand  and  on  order  at  the  beginning  of  period  n, 
prior  to  the  ordering  decision 

y = inventory  on  hand  and  on  order  after  the  ordering  decision  in 
period  n but  prior  to  receiving  demand  in  period  n 
n = number  of  periods  that  the  item  will  be  stocked 
y(n)  = the  appropriate  demand  density  index  in  the  cycle  when  n 
periods  remain 


Mz) 


0 for  z = 0 

1 for  z t 0 


for  y -0 


for  y<0 


V-)  =-  °- 

Def ine 

00 

r.n<y)  - H(y,Y<„»  ♦ - (>Vn)<£> 

S - smallest  integer  that  minimizes  G (•)  . 
n n 


s - smallest  integer  such  that  G (s  ) <K+G  (S). 
n n n n n 

Then 
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and  thus,  ^sn»^n)  Is  an  optimal  policy  when  n periods  remain. 

■ Our  algorithm  calculates  (s  ,S  ),  for  n-A+1,  A + 2 .... 

recursively.  Due  to  the  cyclic  nature  of  the  demand  structure,  an 
optimal  inventory  policy  consists  of  L,  possibly  distinct,  (s,S) 
policies:  that  is,  one  for  each  period  in  the  cycle.  We  terminate 

the  algorithm  when  no  further  change  in  the  policies  occur.  The  just- 
ification for  employing  this  approach  is  mainly  heuristic,  since  we  are 

I 

not  aware  of  any  convergence  theorem  for  this  recursive  process.  In 
fact,  the  computer  program  used  to  implement  the  algorithm  occasionally 
does  not  converge  to  a set  of  policies  but  cycles  between  sets  of 
policies.  An  example  is  given  in  Table  2.1,  which  shows  cycling 

between  94  and  95.  Note,  however,  that  the  difference  in  total  cost 
per  period  between  these  two  inventory  policies  is  negligible. 

\ The  operating  characteristics  of  a specified  inventory  policy  are 

calculated  by  extending  the  methods  of  Wagner  (1969),  who  considered 
the  stationary  environment,  to  our  cyclic  environment.  Wagner  (1969) 
solves  a Markov  chain  to  determine  the  stationary  distribution  of 
inventory  on  hand  and  on  order  after  ordering.  These  probabilities  are 
denoted 

r(y) , for  y « s S, 

where  (s,S)  is  a stationary  inventory  policy  in  use. 

i When  demand  distributions  vary  periodically,  another  dimension  must 

be  added  to  the  Markov  state  space,  namely,  the  period  in  the  cycle.  Thus, 
our  state  space  is  (y,i)  and  represents  the  state  of  having  y units  of 
inventory  on  hand  and  on  order  after  ordering  at  the  beginning  of  the  ltl' 
period  of  the  demand  cycle.  Since  the  policies  used  each  period  may  vary, 
y can  assume  values  from  s^  to  S,  where  s^  is  the  reorder  point  used 
in  the  ith  period  of  the  cycle  and  S is  the  maximum  of  all  L 
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periods  in  the  cycle.  In  order  to  calculate  replenishment  probabilities, 
the  states  (S^i),  for  i = 1,...,L,  are  decomposed  into  states 
(S^(order) , i)  and  S^no  order),!),  for  i * 1,...,L.  The  state 
(S^ (order) , i)  represents  entering  period  i with  inventory  on  hand 
and  on  order  before  ordering  less  than  s^,  so  that  an  order  occurs  in 
period  i.  The  state  (S^no  order),!)  represents  already  having 
units  of  goods  on  hand  and  on  order  before  ordering,  so  that  no  order 
occurs  in  period  i. 

Due  to  the  cyclic  nature  of  the  model,  the  probability  of  entering 
state  (•,.)),  for  Is  j S L,  given  that  the  last  state  is  (*,1),  for 
1 £ i S L,  is  zero  for  j i i + 1,  unless  1 = L and  j * 1.  Thus,  the 
general  structure  of  the  transition  matrix  may  be  presented  by  examining 
only  the  transition  probabilities  from  any  period  i to  period  i + 1, 
for  Is  i S L - 1,  within  the  cycle.  Such  a matrix  is  shown  in  Figure  2.1. 


The  structure  depicted  in  Figure  2.1  assumes  that  s^  < s^  + 

< ^i  + j.  Note  in  Figure  2.1  that  the  structure  of  this  matrix  is  simi- 
lar to  an  upper  triangular  structure,  when  the  last  si  + i ~ rows  are 

ignored.  If  s^  > s(  + , then  the  upper  triangular  structure  would  appear 

only  in  the  first  S - s^  + 1 rows  of  the  matrix. 

Let 

r(Sj(no  order),!)  for  i = 1,...,L 
r (S^ (order) , i)  for  i * 1,...,L 
r(y,i)  for  i =■  1,...,L  and 


y ^ i>1( 

■ 


denote  the  stationary  probabilities  for  the  Markov  chain.  For  noiational 
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rtS^.i)  = r(Sj (order) , 1)  + r(Sj(no  order), i)  for  1 = 

After  solving  the  Markov  chain  for  rC*,*),  we  can  compute  the 
operating  characteristics: 

L S 

Expected  Stockout  Quantity  = E E r(y,i)[  E (Q-y ) (Q)  ] 

1=1  y=s^  Q>y 

L s y X+1 

Expected  Period-End  Inventory  = Z Z r(y,i)[  Z (y-Q)<f>.  (Q)  ] 

i=l  y=si  Q=0 

L 

Replenishment  Frequency  = Z r(S  (order), i) 

i=l  1 

L S y 

Stockout  Probability  =1-1  Z r(y,i)[  E < p A(y-j)] 

i=l  y=s1  j-0 

Expected  Replenishment  Cost  = K(Replenishment  Frequency) 

Expected  Holding  Cost  = h(Expected  Period-End  Inventory) 
Expected  Backlog  Cost  = n (Expected  Stockout  Quantity) 
Expected  Total  Cost  = Expected  Replenishment  Cost  + 

Expected  Stockout  Cost  + 

Expected  Holding  Cost. 

Note  that  each  of  the  L periods  in  the  cycle  contributes  to  the  ex- 
pected value  of  the  operating  characteristic.  Hence,  when  we  refer  to, 
for  example,  expected  total  cost,  we  are  not  referring  to  the  expected 
cost  incurred  in  a specific  period  of  the  cycle,  but  rather  the  average 
value  per  period  for  an  arbitrary  period  in  the  cycle  over  the  infinite 
horizon. 

In  the  stationary  environment  Veinott  and  Wagner  (1965)  prove  that 
if  an  (s,S)  policy  is  optimal,  the  stockout  probability  is  bounded  above 
by  h/(h  + it).  This  upper  bound  is  valid  also  in  the  cyclic  demand 
environment,  as  shown  next. 
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Let  (s^  + C,  + C),  for  1 = represent  an  inventory  policy 

in  a cyclic  demand  environment  of  cycle  length  L.  The  policy 

t"  h 

(s^  + C,  Sj  + C)  is  used  in  the  i period  of  the  cycle.  The  total 
expected  cost  of  such  an  inventory  policy  over  an  entire  cycle  can  be 
expressed  as 

L L D. 

(1)  K 1 r(S  + C (order) , i)  ->  l I1  H(S  + C - y,i)r(S  + C - y>i) 

1=1  1=1  y=0 

where 

S = max  {S^} 


0.  = S - s^,  for  i = 1,...,L. 


We  seek  an  integer  C such  that  total  cost  with  respect  to  the  policy 
(Sj  + C,  Sj  + C),  for  1 = 1,...,L,  is  minimized.  Since  shifting  every 
(s,S)  policy  in  the  cycle  by  C can  not  change  replenishment  frequency, 
minimizing  (1)  with  respect  to  C is  equivalent  to  minimizing 

L D 

Z Z H(S  + C - y , i)r (S  + C - y,i) 
i=l  y=0 

with  respect  to  C. 

It  is  easily  verified  that  H(S  + C - y,i)  is  convex  in  C.  Thus  we 
seek  the  smallest  C such  that 
L D. 

Z i H(S  + C + 1 - y,i)r  (S  + C + I - y.i)  - 
1=1  y-0  l 

Z Z H(S  + C - y,i)r  (S  + C - y.i)  2 0, 
i-1  y-0 

where  r^(*,*)  and  r^(*,‘)  are  the  stationary  distributions  associated  with 
policies  (Sj  + C + 1,  + C + 1),  for  i ■ 1,...,L,  and  (s^  + C,  + C), 

for  i **  1.....L,  respectively.  But  since  shifting  all  the  (s,S)  policies 
In  the  cycle  by  a constant  can  not  affect  the  stationary  distributions 
Tj  (z  + 1,  i)  and  r^  (z,  i)  are  equivalent,  for  z « s^  +C.....S  + C and  for 
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all  i.  Hence,  we  can  use  r'(z,*)  to  denote  r^(z  + 1,  . ) and  r (z,*). 

Velnott  and  Wagner  (1965)  express  H(S  + C + 1 - y,  i) 

- H (S  + C - y,  1)  as  . ...  _ 

(h  + (S  + C - y , i)  -n. 


where 


.A+l  * .A+l,  . 

4>  (z)  * Z (v) 

v-0  1 


Therefore,  we  seek  the  smallest  integer  C such  that 


Z Z1  {(h  + ir)  4^+1(S  + C - y,  i)r  (S  + C - y,  i)  - *r'(S  + C - y,  1)}  2 0 
1*1  y-0  1 

or 


” t 1X1  

Z ZA  4 (S  + C - y,  i)r ' (S  + C - y,  1)  > tt/  (h  + tt). 
i*l  y*0 

The  left  hand  side  of  the  above  inequality  is,  for  the  minimizing  value 
of  C,  the  probability  of  having  0 or  more  units  of  inventory  in  stock. 
Hence,  the  stockout  probability  is  bounded  above  by  1 - it/ (h  + it)  or 
h/(h  + u). 

Since  this  proof  is  valid  for  arbitrary  (s,S)  policies,  it  is 

< t 

valid  for  the  class  of  policies  (s^  + C,  + C),  for  i * 1,...,L, 

f t 

where  (Sj  , S ),  for  i » 1,...,L,  differ  from  the  optimal  policy  by 
C.  Hence,  the  proof  remains  valid  for  optimal  policies. 

2 . 2 Optimal  Control  with  Full  information:  Multi-Item  Systems 

We  let  each  item  in  the  16-item  systems  described  in  Section  1.2.2 
be  controlled  by  using  optimal  replenishment  policies.  Table  2.2  lists 
the  resulting  expected  values  of  average  total  cost  per  period  and  Its 
components  for  the  16- item  negative  binomial  systems  with  variance-to- 
mean  ratio  of  3.  Each  component's  percent  of  total  cost  is  shown  In 


parentheses . 


Table  2.2 
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Average  Costa  Per  Period  for  a Multi-Item 
Negative  Binomial  System  Optimally  Controlled 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

II 

MODEL 

III 

MODEL 

IV 

INVENTORY 

439  (57.8) 

433  (58.0) 

438  (58.2) 

429  (59.7) 

434  (58.6) 

BACKLOG 

109  (14.4) 

106  (14.2) 

106  (14.2) 

102  (14.1) 

104  (14.1) 

REPLENISHMENT 

211  (27.8) 

208  (27.8) 

208  (27.6) 

188  (26.1) 

202  (27.2) 

760  (100.0) 

747  (100.0) 

752  (100.0) 

718  (100.0) 

739  (100.0) 

Observe  that  the  stationary  model's  total  cost  serves  as  an  upper 
bound  for  the  nonstationary  models.  Appendices  A to  E reveal  that  the 
total  cost  of  each  individual  case  in  the  stationary  system  also  provides 
an  upper  limit  for  the  corresponding  cases  in  the  nonstationary  systems. 
Further  experimentation  is  required  to  determine  if  stationary  costs  can 
serve  as  useful  upper  bounds  for  nonstationary  models.  Models  111  and 
IV,  which  have  the  largest  variation  in  demand  means,  have  the  lowest  re- 
plenishment and  total  costs  of  all  the  models  examined. 

The  apportionment  of  expected  average  costs  for  various  classifica- 
tions of  the  items  in  each  multi-item  system  is  shown  in  Table  2.3  which 
gives  the  distributions  by  percentage  of  the  expected  average  total  cost 
per  period.  Aggregate  values  of  other  operating  characteristics  for  the 
various  classifications  are  set  out  in  Table  2.4.  The  classifications  are 
one-way  in  the  sense  that  items  in  the  system  are  grouped  according  to  the 
value  taken  by  a single  input  parameter. 

Table  2.3  reveals  that  the  relative  distributions  of  cost  are  virtu- 
ally the  same  for  the  five  models  that  we  studied.  This  suggests  that  the 
cost  parameters,  the  leadtime,  and  overall  demand  mean  are  the  dominant 
factors  in  determining  the  distributions  of  system  cost,  whereas  the  na- 
ture of  the  nonstationarlty  has  little  inf  u nee  on  the  cost  distributions. 
With  the  exception  of  the  replenishment  fre  :'iPncieS  of  models  Hi  and  IV, 


Table 
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Apportionment  of  Aggregate  Average  Costs  Per  Period  for  a Multi-Item  System  Under 
Optimal  Control  with  Full  Information,  for  Various  Classifications 


Apportionment  of  Aggregate  Average  Costs  Per  Period  for  a Multi-Item  System  Under 


TOTAL  | 100.0  38.7  61.3  45.5  54.5  47.2  52.8  ; 42.0  57.9 


Apportionment  of  Aggregate  Average  Costs  Per  Period  for  a Multi-Item  System  Under 
Optimal  Control  with  Full  Information,  for  Various  Classifications 
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Table  2.4a 


Operating  Characteristics  of  a Multi-Item  System  Under 
Optimal  Control  with  Full  Information,  for  Various  Classifications 


3 


i 

v£> 


REPLENISHMENT  FREQUENCY  | .284  .270  .299  .326  .243  .290  .279  .232  .337 


Table  2.4c 
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the  aggregate  values  of  operating  characteristics  in  Table  2.4  of  the 
nonstationary  systems  are  similar  to  those  of  the  stationary  model. 

The  sensitivity  of  operating  characteristics  and  policy  parameters 
is  summarized  in  Table  2.5,  which  has  been  constructed  by  observing  any 
monotonicity  in  the  relation  for  a multi-item  system.  With  the  exception 
of  the  stationary  model’s  decreasing  relationship  of  backlog  frequency 
with  leadtime,  the  relationships  of  operating  characteristics  to  input 
parameters  are  identical  for  the  stationary  system  and  all  non-station- 
ary  systems,  suggesting  that  the  nature  of  the  nonstationarity  does  not 
play  a major  role  in  the  sensitivity  of  operating  characteristics  to  in- 
put parameters. 

Policy  parameters  s and  S also  follow  identical  relationships 
for  those  systems  studied;  however,  the  difference  S-s,  which  is  denoted  by 

D,  does  not  generally  maintain  the  same  relationships  for  each  system.  For 
example,  while  D increases  with  respect  to  leadtime  in  the  stationary 
system,  no  such  simple  relationship  exists  in  any  of  the  nonstationary 
systems.  The  only  monotonic  relationship  for  D that  all  systems  have 
in  common  is  that  with  respect  to  replenishment  cost. 

Finally,  we  explore  one  ofher  point,  an  approximation  well-known 
to  writers  on  inventory  control,  which  is  potentially  of  great  practical 
use.  As  proved  in  Section  2.1,  for  an  optimally-controlled  multi-period 
inventory  model  with  linear  costs,  the  probability  of  a stockout  occurlng 
is  no  larger  than  h/(h  + rr) . Inspection  of  Tables  2.4  shows  that  the 
backlog  frequencies  are  close  to  their  upper  bounds  for  every  system.  For 
example,  the  8 items  having  it  - 4 in  the  model  I system  have  a combined 
relative  backlog  frequency  of  .192,  compared  to  the  upper  bound  of  .200. 
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Table  2.5a 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics 
of  an  Optimally  Controlled  Inventory  System  When 
Demand  Has  a Negative  Binomial  Distribution 

INPUT  PARAMETERS 

C(OUT)/C(IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 

STATIONARY  MODEL 


OPERATING  CHARACTERISTIC 

E(COST)  + + 

E( PERIOD-END  INVENTORY)  f t 

E( BACKLOG  COST)  1 + 

E( BACKLOG  FREQUENCY)  + + 

PROPORTION  DEMAND  BACKLOGGED  1 + 

E( REPLENISHMENT  COST)  + t 

E (REPLENISHMENT  FREQUENCY)  + + 

POLICY  PARAMETERS 

D + t 

s t + 

S t t 

MODEL  I 

OPERATING  CHARACTERISTIC 

E(COST)  t t 

E(PERIOD-END  INVENTORY)  f t 

E( BACKLOG  COST)  1 t 

E( BACKLOG  FREQUENCY)  + + 

PROPORTION  DEMAND  BACKLOGGED  + t 

E( REPLENISHMENT  COST)  f t 

E (REPLENISHMENT  FREQUENCY)  t + 

POLICY  PARAMETERS 

D + + 

s t 1 

S t t 


t indicates  relation  is  monotone  increasing 
1 indicates  relation  is  monotone  decreasing 
^ indicates  relation  is  approximately  constant 
? indicates  no  conclusion  can  be  drawn  from  the  dala 
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Table  2.5b 

Sensitivity  of  Input  Parameters  of  Operating  Characteristics 
of  an  Optimally  Controlled  Inventory  System  When 
Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT)/C(IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 


STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 

E(COST) 

t 

F 

F 

F 

E(PERIOD-END  INVENTORY) 

f 

F 

F 

F 

E( BACKLOG  COST) 

F 

F 

F 

F 

E( BACKLOG  FREQUENCY) 

F 

F 

F 

F 

PROPORTION  DEMAND  BACKLOGGED 

F 

F 

F 

F 

E (REPLENISHMENT  COST) 

F 

F 

F 

F 

E (REPLENISHMENT  FREQUENCY) 

F 

F 

F 

F 

POLICY  PARAMETERS 

D 

F 

F 

F 

F 

s 

F 

F 

F 

t 

s 

F 

F 

F 

F 

MODEL  II 

OPERATING  CHARACTERISTIC 

E(COST) 

t 

F 

F 

F 

E (PERIOD- END  INVENTORY) 

t 

F 

F 

F 

E ( BACKLOG  COST) 

F 

F 

F 

F 

E (BACKLOG  FREQUENCY) 

F 

F 

F 

F 

PROPORTION  DEMAND  BACKLOGGED 

F 

F 

F 

F 

E (REPLENISHMENT  COST) 

F 

F 

F 

F 

E (REPLENISHMENT  FREQUENCY) 

t 

F 

F 

F 

POLICY  PARAMETERS 

D 

7 

F 

7 

? 

8 

F 

F 

F 

F 

S 

F 

F 

F 

F 

F indicates  relation  is  monotone  increasing 
1 indicates  relation  is  monotone  decreasing 
'v  indicates  relation  is  approximately  constant 
? indicates  no  conclusion  can  be  drawn  from  the  data 
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Table  2.5c 


t 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics 
of  an  Optimally  Controlled  Inventory  System  When 
Demand  Has  a Negative  Binomial  Distribution 


C (OUT) /C (IN ) 


STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 

E(COST)  t 

E(PER10D-END  INVENTORY)  t 

E( BACKLOG  COST)  + 

E (BACKLOG  FREQUENCY)  + 

PROPORTION  DEMAND  BACKLOGGED  I 

E (REPLENISHMENT  COST)  t 

E (REPLENISHMENT  FREQUENCY)  t 

POLICY  PARAMETERS 

D 1 

s t 

S t 

MODEL  III 

OPERATING  CHARACTERISTIC 

E(COST)  t 

E (PERIOD-END  INVENTORY)  t 

E( BACKLOG  COST)  1 

E( BACKLOG  FREQUENCY)  + 

PROPORTION  DEMAND  BACKLOCCED  1 

E( REPLENISHMENT  COST)  1 

E( REPLENISHMENT  FREQUENCY)  t 

POLICY  PARAMETERS 

D ? 

s t 

S 1 


INPUT  PARAMETERS 
C(FIX)/C(IN)  LEADTIME 


+ indicates 
1 indicates 
v indicates 
? indicates 


relation  is  monotone  increasing 
relation  is  monotone  decreasing 
relation  is  approximately  constant 
no  conclusion  can  be  drawn  from  the  data 


MEAN  DEMAND 


t 

t 

t 

1 

+ 

t 

t 


t 

1 

t 


t 

t 

t 

t 

1 

t 

t 


? 

t 

t 
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Sensitivity  of  Input  Parameters  of  Operating  Characteristics 
of  an  Optimally  Controlled  Inventory  System  When 
Demand  Has  a Negative  Binomial  Distribution 


C(0UT)/C(1N) 


STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 

E(COST)  t 

E (PERIOD- INVENTORY)  t 

E( BACKLOG  COST)  + 

E( BACKLOG  FREQUENCY)  1 

PROPORTION  DEMAND  BACKLOGGED  + 

E (REPLENISHMENT  COST)  t 

E (REPLENISHMENT  FREQUENCY)  t 

POLICY  PARAMETERS 

D 1 

s + 

S + 

MODEL  IV 

OPERATING  CHARACTERISTIC 

E(COST)  f 

E (PERIOD-END  INVENTORY)  f 

E( BACKLOG  FREQUENCY)  i 

PROPORTION  DEMAND  BACKLOGGED  1 

E (REPLENISHMENT  COST)  f 

E (REPLENISHMENT  FREQUENCY)  f 

POLICY  PARAMETERS 

D ? 

8 + 

S t 


INPUT  PARAMETERS 
C(FIX)/C(IN)  LEADTIME 


1 indicates  relation  is  monotone  increasing 
* indicates  relation  is  monotone  decreasing 
^ indicates  relation  is  approximately  constant 
? Indicates  no  conclusion  can  be  drawn  from  the  data 


MEAN  DEMAND 


t 

t 

t 

t 

i 

t 

t 


t 

t 

t 


t 

t 

1 

t 

I 

I 

t 

t 

t 

t 


I 
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2 . 3 Optimal  Control  with  Full  Informa t Ion : Single-I  tern  Syst rms 

Summary  information  about  the  performance  of  single-item  systems  under 
optimal  control  with  full  information  is  tabulated  in  Appendices  A to  E. 

Two  single-item  cases  are  selected  for  the  collection  of  more  detailed  in- 
formation. Operating  costs  are  examined  for  each  period  of  the  12  month 
cycle.  For  convenient  reference,  the  cases  are  labelled  A and  B,  and  the 
input  parameters  are  set  out  in  Table  2.6.  These  two  cases  are  chosen  be- 
cause they  best  illustrate  the  effect  of  the  nonstationary  demand  structure 
on  the  operating  characteristics.  To  facilitate  the  analysis.  Figure  2.2 

displays  the  A + 1 fold  convoluted  demand  means  for  both  cases.  Note  that 

_( A+l ) 

the  convoluted  means  differ  from  only  by  the  factor  1/(A  + 1)  and 

a . a v.  (*+1) 

are  denoted  by 

Table  2.6 

2 Single-Item  Cases 


Case 

! Demand 
i Structure 

MEAT 

! 

LEADTIME  j 

| 1 

C (FIX) /C (IN) 

C (OUT) /C ( IN) 

A 

Model  I 

16 

2 

32 

99 

B 

Model  IV 

16 

2 

32 

99 

Figure  2.3  shows  the  monthly  cost  contributions  to  average  total  cost 
per  period  and  its  cost  components  for  case  A.  Replenishment  costs  for 
month  I are  the  result  of  orders  made  at  the  beginning  of  month  i. 

Since  an  order  decision  at  the  beginning  of  month  l can  not  affect  the 
system  until  the  order  arrives  A periods  later,  we  associate  the  backlog 
and  holding  costs  incurred  at  the  end  of  period  l + A with  period  i. 

Thus,  backlog  and  holding  costs  graphed  in  Figure  2.3  actually  occur  A 
periods  later.  Table  2.7  lists  case  A's  optimal  (s^.S^)  policy,  for 
1 - 1 12. 


I 
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mean 
Of MONO 
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°FRI0D  IN  CYCLE 


a.  Case  A 


“1 1 1 f 1 1 1 1 1 1 1 1 

1.00  2.00  3.00  4 00  5 00  6 00  7 00  8.00  9.QQ  10  0 11  0 12.0 

PERIOD  IN  CYCt  F 
b.  Case  B 


Figure  2.2  Convoluted  Mean  Demands  for  Case  A and  Case  P 
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The  most  striking  property  of  Figure  2.3  is  the  presence  of  two 
large  spikes  in  each  of  the  cost  components.  This  phenomenon  is  best 
explained  by  examining  each  period's  order  probability,  which  is  propor- 
tional to  its  replenishment  costs,  and  the  optimal  policy. 

Table  2.7  reveals  that  s^  is  relatively  stable  for  the  first 
seven  periods  of  the  cycle.  Similarly,  is  relatively  constant  for 

the  fitst  five  periods  of  demand  but  begins  to  decrease  in  periods  6 and 
7.  The  policy  employed  in  period  8 deviates  sharply  from  this  pattern. 

The  policy  parameters  s and  S increase  greatly  over  previous  values. 

o o 

This  is  caused  by  the  value  of  = 25.6  , which  is  twice  the  value  of 

(3) 

the  first  nine  demand  means  in  the  cycle.  Thus,  PQ  is  the  first  convo- 

o 

luted  demand  mean  in  the  cycle  to  reflect  the  presence  of  larger  demands. 
In  anticipation  of  the  increased  demand,  Sg  greatly  exceeds  the  first 
seven  values  of  s.  in  the  cycle,  so  that  the  probability  of  ordering 

(3) 

in  period  8 increases  sharply.  As  p ^ increases,  so  does  sJt  both 

reaching  maximum  values  in  period  10.  The  parameter  also  increases 

(3) 

roughly  with  , but  the  maximum  value  for  occurs  in  period  9, 

one  period  before  the  largest  convoluted  mean  is  encountered.  Hence,  the 

movement  of  optimal  policies  for  case  A can  be  analyzed  by  examining  the 

behavior  of  the  convoluted  means,  shown  in  Figure  2.2. 

Since  s„  is  several  units  greater  than  S7,  there  is  a high  proba- 

bility  of  ordering  in  period  8.  The  increase  of  S over  past  values 

o 

results  in  an  increase  in  expected  inventory  costs  and  a corresponding 
decline  in  expected  penalty  costs.  If  we  assume  that  an  order  has  taken 
place  in  period  8,  then  we  enter  period  9 with  an  expected  value  of 
108  - 12.8  • 95.2  units  of  inventory  on  hand  and  on  order,  which  greatly 
exceeds  s^  - 86.  Thus,  the  order  probability  in  period  9 decreases. 


causing  a decrease  in  expected  holding  cost  and  an  Increase  in  expected 


TOTAL  COST 


Components  for  Case  A Under  Optimal  Control  With  Full  Information 
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penalty  cost.  Assuming  that  no  order  is  placed  in  period  9,  we  then 
enter  period  10  with  an  expected  95.2  - 12.8  = 82.4  units  of  inventory 
on  hand  and  on  order,  which  is  significantly  less  than  Sj(  = 105.  Hence, 
we  again  see  an  increase  in  order  probability  that  increases  expected 
holding  costs  and  reduces  expected  backlog  costs.  After  period  10's 
demand  the  expected  on  hand  and  on  order  level  is,  assuming  an  order 
occurs  in  period  10,  S1Q  - 25.6  - 102.4. 

It  is  important  to  note  that  in  explaining  the  structure  of  the 

cost  components  in  the  demand  cycle,  we  have  assumed  a deterministic 
pattern  in  the  order  policy.  Nevertheless,  this  deterministic  approach 
approximates  the  actual  probabilistic  model  well  enough  to  explain  the 
presence  of  spikes  in  the  cost  components.  If  we  were  to  continue  the 
analysis,  we  would  observe  a small  probability  of  ordering  in  period  11 
followed  by  a larger  probability  of  ordering  in  period  12.  After  period 
12  the  (s,S)  policies  and  mean  period  demands  are  stable.  This  stability 
is  reflected  in  the  cost  components  which  are  approximately  constant  in 
periods  1 through  6. 

The  underlying  mean  demand  structure  of  case  B,  shown  in  Figure 
1.1. e,  exhibits  considerably  more  gradual  change  than  is  observed  in 
case  A.  This  slower  movement  results  in  a gradual  change  in  the  opti- 
mal policies  and  their  cost  components.  Table  2.8  displays  the  optimal 
policy  for  case  B.  Note  that  s.  and  attain  their  lowest  values 

in  periods  2 and  3 and  highest  values  in  periods  8 and  9.  The  convo- 
(3) 

luted  mean,  . plotted  in  Figure  2.2.b,  attains  its  minimum  and 

maximum  values  in  the  same  periods,  respectively.  The  policy  para- 
meters Sj  and  follow  a sinusoidal  pattern,  similar  to  that  of 


The  12  period  cost  components,  displayed  in  Figure  2.4,  also  follow 

(3) 

the  sinusoidal  pattern  of  . The  cost  components  do  not  reflect  the 
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shape  of  as  well  as  s^  and  do,  but  they  nevertheless  echo 

the  basic  shape  of  the  convoluted  means.  Backlog  costs  deviate  the 
most  from  this  pattern;  even  for  this  component,  however,  its  minimum 
and  maximum  values  occur  near  the  corresponding  extreme  points  of  . 

Unlike  case  A,  case  B does  not  display  significant  spiking  in  any 
of  its  operating  characteristics.  Case  B's  mean  demand  structure  con- 
tains gradual  movements,  in  contrast  to  the  step  change  of  case  A.  Thus, 
while  case  B's  reorder  policy  can  adjust  slowly  to  increasing  demand, 
case  A must  adjust  immediately  to  the  first  indication  of  higher  mean 
demands,  producing  spikes  in  the  cost  components. 

2 . 4 Approximately  Optimal  Control  with  Full  Information:  Multi-item 

Systems 

The  replenishment  policy  parameter  values  employed  are  calculated 
using  the  modified  version  of  the  Ehrhardt  power  approximation  (1976) 
discussed  in  Section  1.1.2.  Although  the  control  is  described  as 
utilizing  full  information,  actually  only  the  means  and  variances  of 
the  demand  distributions  need  to  be  known  in  order  to  compute  an  inven- 
tory policy  of  the  (s,S)  type. 

The  components  of  aggregate  total  coat  per  period  for  the  various  multi 
item  systems  are  recorded  in  Table  2.9,  and  Table  2.10  contains  a com- 
parison of  the  systems  controlled  optimally  and  approximately  optimally. 

Note  that  for  all  systems  the  power  approximation  results  in  only 
small  increases  of  total  cost  above  the  optimal  total  cost.  The  largest 
increases  occur  under  those  models  with  the  greatest  variation  in  demand 
means:  model  III  and  model  IV.  For  these  models  the  increases  in  total 
cost  are  5.8%  and  2.6%,  respectively.  Table  2.9  also  reveals  that  the 
inventory  costs  of  each  nonstationary  model  exceeds  that  of  the  station- 
ary model,  while  nonstationary  backlog  and  replenishment  costs  are  lower 
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than  those  of  the  stationary  model.  These  differences  are  most  notice- 
able when  contrasting  models  III  and  IV  with  the  stationary  model. 


Tabic  2.9 

Average  Coata  Per  Period  for  a Multi-Item 
Negative  Binomial  System  Controlled  Approximately  Optimally 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

11 

MODEL 

III 

MODEL 

IV 

INVENTORY 

440  (57.8) 

445  (58.8) 

442  (58.4) 

466  (61.3) 

461  (60.7) 

BACKLOG 

109  (14.3) 

105  (13.8) 

106  (14.0) 

95  (12.5) 

97  (12.8) 

REPLENISHMENT 

212  (27.8) 

207  (27.3) 

209  (27.6) 

199  (26.2) 

201  (26.5) 

TOTAL 

761  (100.0) 

757  (100.0) 

758  (100.0) 

760  (100.0) 

759  (100.0) 

Note:  Numbers  In  parentheses  are  percent  of  total. 


Table  2.10 

Average  Costa  Per  Period  for  a Multi-Item 
Negative  Binomial  System 

Comparison  of  Optimal  and  Approximately  Optimal  Control 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

II 

MODEL 

III 

MODEL 

IV 

INVENTORY 

1 (.1) 

12  (2.9) 

5 (1.1) 

37  (8.7) 

27  (6.3) 

BACKLOG 

0 (-.2) 

-1  (-1.0) 

0 (-.2) 

-7  (-6.5) 

-7  (-7.0) 

REPLENISHMENT 

1 (.3) 

-1  (-.5) 

1 (.5) 

11  (6.0) 

0 (-.2) 

TOTAL 

1 (-U 

10  (1.4) 

6 (.8) 

42  (5.8) 

19  (2.6) 

Note:  Table  2.10  shows  the  absolute  Increase  or  decrease  In  the  coat  components  of 
Tuble  2.9  over  those  In  Table  2.2  with  percent  changes  in  parentheses. 


Table  2.10  indicates  that  the  nonstationary  systems  maintain  higher 
levels  of  inventory  than  that  found  in  the  corresponding  optimal  systems 
which  result  in  lower  backlog  costs  under  approximately  optimal  control. 
This  is  most  apparent  in  models  III  and  IV. 

The  apportionment  of  aggregate  total  cost  per  period  for  various 
classifications  of  the  items  in  each  system  is  shown  in  Table  2.11,  which 
can  be  compared  with  Table  2.3.  Observe  that  all  systems  follow  approxi- 
mately the  same  apportionment  of  costs.  Model  Ill's  and  model  IV's 
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Apportionment  of  Aggregate  Average  Costs  Per  Period  for  a Multi-Item  System  Under 
Approximately  Optimal  Control  with  Full  Information,  for  Various  Classifications 


Percentage  Excess  of  Costs  Per  Period  of  a Multi-Item  System  Under 
Approximately  Optimal  Control  Over  the  Costs  of  the  Same  System  Under  Optimal  Control 


Percentage  Excess  of  Costs  Per  Period  of  a Multi-Item  System  Under 


57 


TOTAL 
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apportionment  of  inventory  costs,  however,  are  slightly  greater  than 
the  other  systems.  Detailed  comparison  of  the  systems  under  optimal 
and  approximately  optimal  control,  expanding  on  Table  2.10,  is  in 
Table  2.12,  showing  the  percentage  increases  per  period  in  each  com- 
ponent of  cost. 

Tables  2.12a  and  2.12b  show  that  under  low  unit  penalty  costs, 
the  stationary  model,  model  I,  and  model  II  incur  lower  than  optimal 
inventory  costs  at  the  expense  of  higher  backlog  costs.  This  strategy 
is  reversed  when  the  unit  penalty  cost  is  99.  Under  the  higher  penalty 
cost,  these  systems  carry  excess  inventory  to  lessen  the  incidence  of 
backlogs.  Models  III  and  IV,  however,  carry  excess  inventory  regard- 
less of  the  value  of  unit  penalty  cost,  although  small  penalty  costs 
result  in  smaller  excess  inventory  levels  than  do  larger  penalty  costs. 

Table  2.13  gives  the  values  of  other  operating  characteristics 
for  the  systems  under  approximately  optimal  control.  The  values  of 
these  operating  characteristics  for  model  III  and  model  IV  are  different 
from  those  of  the  remaining  systems.  Model  III  and  IV’s  high  Inventory 
levels  are  distributed  differently  from  the  other  systems,  and  their 
backlog  frequencies  are  less  than  those  found  In  the  stationary  model, 
model  I,  and  model  II.  These  differences  are  also  apparent  in  Table 
2.14,  which  gives  the  percentage  changes  in  the  values  of  the  operating 
characteristics  in  Table  2.13  from  the  corresponding  ones  in  Table  2.4. 

Table  2.15  illustrates  the  sensitivity  of  operating  characteris- 
tics and  policy  parameters  to  input  parameters.  With  the  exception  of 
backlog  frequency  with  respect  to  leadtime,  the  operating  character- 
istics of  all  four  nonstationary  systems  follow  Identical  relationships 
over  their  input  parameters.  The  stationary  model,  model  I and  model 


Table  2.13b 


62 


Table  2.13c 
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perating  Characteristics  of  a Multi-Item  System  Under  Approximately 
Optimal  Control  with  Full  Information,  for  Various  Classifications 


REPLENISHMENT  FREQUENCY  .274  .278  .269  .309  .238  .276  .272 


Percentage  Excess  of  Operating  Characteristics  of  a Multi-Item  System  Under  Approximately 
Optimal  Control  Over  the  Same  Characteristics  for  the  System  Under  Optimal  Control 
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Percentage  Excess  of  Operating  Characteristics  of  a Multi-Item  System  Under  Approximately 
Optimal  Control  Over  the  Same  Characteristics  for  the  System  Under  Optimal  Control 
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Percentage  Excess  of  Operating  Characteristics  of  a Multi-Item  System  Under  Approximately 
Optimal  Control  Over  the  Same  Characteristics  for  the  System  Under  Optimal  Control 


6 


REPLENISHMENT  FREQUENCY  7.1  12.8  1.8  10.0  3.5  7.4  6.7  3.6  9.4 


Percentage  Excess  of  Operating  Characteristics  of  a Multi-Item  System  Under  Approximately 
Optimal  Control  Over  the  Same  Characteristics  for  the  System  Under  Optimal  Control 
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IV  have  increasing  backlog  frequencies  with  respect  to  leadtime,  whereas 
models  II  and  III  have  decreasing  backlog  frequencies.  Note  also  that 
unlike  the  stationary  model,  the  nonstationary  models  do  not  follow  con- 
stant relationships  with  respect  to  unit  penalty  costs  for  replenish- 
ment costs  and  replenishment  frequency. 

Policy  parameters  s and  S follow  identical  relationships  in 
all  systems  under  approximately  optimal  control.  These  relationships 
also  are  in  agreement  with  those  observed  under  optimal  control. 

Policy  parameter  D,  under  approximately  optimal  control,  follows  an 
increasing  relationship  with  respect  to  mean  demand  in  all  systems, 
except  in  model  III  for  those  items  with  a replenishment  cost  of  32 
and  a leadtime  of  2.  Further  examination  of  model  III  reveals  that 
the  largest  values  for  p in  all  the  nonstationary  systems  occur 

in  that  model.  The  power  approximation  first  computes  D using  equa- 
tion (1)  in  Section  1.1.2;  D/p^**^  is  a measure  of  order  frequency, 
and  if  this  quantity  is  less  than  or  equal  to  1.5,  the  algorithm 
attempts  to  find  a smaller  value  for  D that  will  increase  order 
frequency  [see  Ehrhardt  (1976)].  Under  small  leadtimes  and  unit 
replenishment  costs,  model  Ill's  demand  structure  occasionally  resulted 
in  small  values  for  D that  were  not  monotone  with  respect  to  mean 
demand  for  large 

Finally,  we  examine  the  ratio  1/(1  + *)  as  an  approximation  to  the 
probability  of  backlog  in  any  period  for  each  item  (recall  h « 1) . The 
entries  in  Table  2.16  have  been  extracted  from  Table  2.13.  The  approxi- 
mation is  a reasonable  one  for  the  stationary  model  and  models  1 and  II. 
For  models  III  and  IV,  however,  the  actual  backlog  probability  is 
slightly  overstated  by  1/(1  + it). 


Table  2.15a 
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Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of 
an  Approximately  Optimally  Controlled  Inventory  System 
When  Demand  Has  a Negative  Binomial  Distribution 

INPUT  PARAMETERS 

C(0UT)/C(IN)  C(FIX)/C(1N)  LEADTIME  MEAN  DEMAND 

STATIONARY  MODEL 


OPERATING  CHARACTERISTIC 

E(COST)  t + 
E(PERIOD-END  INVENTORY)  t 1 
E( BACKLOG  COST)  t 
E( BACKLOG  FREQUENCY)  + t 
PROPORTION  DEMAND  BACKLOGGED  + t 
E( REPLENISHMENT  COST)  -v  t 
E (REPLENISHMENT  FREQUENCY)  -v  + 


t 

+ 

t 

t 

i 

* 

+ 


POLICY  PARAMETERS 
D 


% 


t 


s 


t 


S 


+ 


MODEL  I 

OPERATING  CHARACTERISTIC 


E(COST)  t 
E (PERIOD-END  INVENTORY)  t 
E( BACKLOG  COST)  + 
E( BACKLOG  FREQUENCY)  + 


t 

t 

t 


i 

+ 

t 

+ 


PROPORTION  DEMAND  BACKLOGGED  + + 

E( REPLENISHMENT  COST)  + t 

E (REPLENISHMENT  FREQUENCY)  + + 

POLICY  PARAMETERS 

D 'v  t 

s t + 

S t 1 


f 

+ 

+ 

7 

t 

f 


i indicates  relation  is  monotone  increasing 
T indicates  relation  is  monotone  decreasing 
'v  indicates  relation  is  approximately  constant 
? indicates  no  conclusion  can  be  drawn  from  the  data 


t 

t 

t 

+• 

+ 

t 

t 


t 

t 

t 


t 

+ 

t 

+ 

+ 

t 

t 

t 

t 

i 
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Table  2.15b 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of 
an  Approximately  Optimally  Controlled  Inventory  System 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT) /C(IN)  C(FIX) /C(IN)  LEADTIME 

STATIONARY  MODEL 


OPERATING  CHARACTERISTIC 


E(COST)  t 4 

E( PERIOD-END  INVENTORY)  t 4 

E (BACKLOG  COST)  4 4 

E (BACKLOG  FREQUENCY)  4 t 

PROPORTION  DEMAND  BACKLOGGED  4 4 

E( REPLENISHMENT  COST)  'v  4 

E( REPLENISHMENT  FREQUENCY)  % 4 

POLICY  PARAMETERS 

D % 4 

s t 1 

S it 

MODEL  II 

OPERATING  CHARACTERISTIC 

E(COST)  4 4 

E( PERIOD-END  INVENTORY)  t 4 

E( BACKLOG  COST)  4 t 

E( BACKLOG  FREQUENCY)  1 4 

PROPORTION  DEMAND  BACKLOGGED  4 t 

E (REPLENISHMENT  COST)  4 4 

E( REPLENISHMENT  FREQUENCY)  4 4 

POLICY  PARAMETERS 

D -vi 


4 

4 

4 

t 

t 

4 

4 

4 

t 

4 


4 

4 

4 

4 

t 

4 

4 

? 


8 


4 


4 


t 

4 

'V 

? 


S 

indicates 

indicates 

indicates 

indicates 


4 4 

relation  is  monotone  increasing 
relation  is  monotone  decreasing 
relation  is  approximately  constant 
no  conclusion  can  be  drawn  from  the  data 


4 


MEAN  DEMAND 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 
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Table  2.15c 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of 
an  Approximately  Optimally  Controlled  Inventory  System 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C (OUT) /C (IN)  C(FIX) /C ( IN)  LEADTIME 


STATIONARY  MU DEL 


OPERATING  CHARACTERISTIC 


E(COST)  4 t 

E (PERIOD-END  INVENTORY)  4 4 

E( BACKLOG  COST)  4 4 

E( BACKLOG  FREQUENCY)  4 4 

PROPORTION  DEMAND  BACKLOGGED  I 4 

E (REPLENISHMENT  COST)  % 4 

E( REPLENISHMENT  FREQUENCY)  'v-  4 

POLICY  PARAMETERS 

D % 4 

s 4 4 

S 4 4 

MODEL  III 

OPERATING  CHARACTERISTIC 

E(COST)  4 t 

E (PERIOD-END  INVENTORY)  4 4 

E ( BACKLOG  COST)  I t 

E( BACKLOG  FREQUENCY)  + 4 

PROPORTION  DIMAND  BACKLOGGET  4 t 

E( REPLENISHMENT  COST)  1 4 

E( REPLENISHMENT  FREQUENCY)  4 4 

POLICY  PARAMETERS 


4 

4 

4 

4 

4 

4 

4 

t 

4 

t 


4 

t 

4 

4 

4 

4 

4 


7 
t 
4. 

4 indicates  relation  is  monotone  increasing 
4 indicates  relation  is  monotone  decreasing 
'v  indicates  relation  is  approximately  constant 
? indicates  no  conclusion  can  be  drawn  from  the  data 
* indicates  an  explanation  is  given  in  the  text 


MEAN  DEMAND 


t 

t 

t 

4 

4 

t 

t 


4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


* 

4 

4 


* 
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Table  2.15d 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of 
an  Approximately  Optimally  Controlled  Inventory  System 
When  Demand  Has  a Negative  Binomial  Distirbution 


INPUT  PARAMETERS 

C(OUT)/C(IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 

STAT IONARY  MODEL 


OPERATING  CHARACTERISTIC 

E(COST)  4 1 

E (PERIOD- END  INVENTORY)  4 4 

E (BACKLOG  COST)  4 4 

E (BACKLOG  FREQUENCY)  1 4 

PROPORTION  DEMAND  BACKLOGGED  1 4 

E (REPLENISHMENT  COST)  'v  4 

E (REPLENISHMENT  FREQUENCY)  -v  1 

POLICY  PARAMETERS 

D 'v*  4 

s 4 4 

S 4 4 

MODEL  IV 

OPERATING  CHARACTERISTIC 

E(COST)  4 4 

E (PERIOD- END  INVENTORY)  4 t 

E( BACKLOG  COST)  4 t 

E( BACKLOG  FREQUENCY)  4 t 

PROPORTION  DEMAND  BACKLOGGED  4 4 

E (REPLENISHMENT  COST)  4 4 

E (REPLENISHMENT  FREQUENCY)  4 4 

POLICY  PARAMETERS 

D n,  t 

s +4 

S it 


t indicates 
4 indicates 
'v  indicates 
7 indicates 


relation  is  monotone  Increasing 
relation  is  monotone  decreasing 
relation  is  approximately  constant 
no  conclusion  can  be  drawn  from  the  data 


4 

t 

+ 

4 

t 

4 

4 


+ 

+ 

t 


t 

4 

t 

t 

+ 

4 

4 

? 

t 

t 


t 

4 

t 

4 

4 

t 

4 


+ 

t 

t 


t 

4 

t 

4 

4 

t 

4 


t 

4 

4 


Table  2.16 


Backlog  Frequencies 
With  Negative 

for  Multi-Item  Systems 
Binomial  Demands 

C(0UT)/C(IN) 

(-»>. 

4 

99 

l/U-Hr) 

.2 

.01 

SYSTEM: 

Stationary  Model 

.200 

.009 

Model  I 

.195 

.009 

Model  II 

.197 

.009 

Model  III 

.178 

.008 

Model  IV 

.181 

.008 

2 . 5 Approximately  Optimal  Control  with  Full  Information: 

Single-Item  Systems 

The  approximate  policies  of  case  A,  displayed  in  Table  2.17,  follow 
the  same  pattern  observed  in  the  optimal  policy.  Case  A's  power  approxi- 
mation policy  maintains  a constant  policy  for  the  first  seven  periods  of 
demand.  Both  and  increase  in  periods  8 through  10  and  decrease 

in  periods  11  and  12,  following  the  similar  pattern  of  the  optimal 
policy.  Comparing  Table  2.17  with  Table  2.7  reveals  that  although  opti- 
mal and  app-oximately  optimal  policies  follow  the  same  basic  patterns, 
the  two  policies'  values  for  may  differ  a great  deal  in  some  periods 

of  the  cycle.  The  total  percent  excess  cost  for  case  A,  however,  is  only 
1.86%,  suggesting  that  acceptable  inventory  policies  may  differ  a good 
deal  from  optimal  policies. 

Figure  2.5  plots  the  cost  components  for  Case  A.  Holding  and  replen- 
ishment costs  govern  the  general  shape  of  the  total  cost  curve.  Unlike 
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Table  2.17 

Approximately  Optimal  Policy  for  Case  A 


Period 
Cycle  : 

in 

i 

s. 

1 

Si 

1 

57 

(57) 

85 

(82) 

2 

57 

(57) 

85 

(83) 

3 

57 

(57) 

85 

(82) 

4 

57 

(57) 

85 

(82) 

5 

57 

(57) 

85 

(84) 

6 

57 

(59) 

85 

(77) 

7 

57 

(55) 

85 

(67) 

8 

72 

(72) 

104 

(108) 

9 

87 

(86) 

122 

(130) 

10 

102 

(105) 

140 

(128) 

11 

87 

(89) 

122 

(112) 

12 

72 

(73) 

104 

(95) 

Note:  The 

optima  1 

policy 

is  given 

in  parentheses. 

t 


0T3L  COST 


Components  for  Case  A Under  Approximately  Optimal  Control  With  Full  Information 
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case  A under  optimal  control,  peaks  in  the  replenishment  cost  function 
do  not  always  produce  peaks  in  backlog  cost.  This  is  most  noticeable 
in  period  8 where  a large  increase  in  reorder  probability  has  little 
effect  on  the  corresponding  backlog  component. 

Table  2.18  lists  the  approximate  policy  for  case  B.  As  under  op- 
timal control,  the  approximate  policy  closely  follows  the  basic  move- 
ment of  the  convoluted  means.  The  differences  between  in  the  opti- 

mal and  power  approximations  policies  are  significant  in  some  periods; 
however,  the  power  approximation  results  in  excess  total  costs  of  only 

1.32%.  The  cost  components  of  case  B,  illustrated  in  Figure  2.6  reflect 

(3) 

the  movement  of  in  the  same  manner  as  the  optimal  cost  components. 

For  each  period  in  the  cycle,  the  total  cost  function  under  approxi- 
mately optimal  control  does  not  deviate  greatly  from  the  total  cost 
function  under  optima]  control,  shown  in  Figure  2.4. 

2 . 6 The  Effects  of  Demand  Structure  Misspecif ication  on  the  Power 
Approximation 

In  Section  2.4  and  2.5  we  assumed  knowledge  of  the  demand  mean  and 
variance  for  each  period  in  the  12  month  cycle.  The  decision  maker, 
however,  does  not  always  have  complete  knowledge  of  the  underlying  de- 
mand structure.  Section  2.6  investigates  the  effects  of  demand  struct- 
ure misspecif ication  on  the  power  approximation. 

We  first  examine  the  situation  in  which  the  decision  maker  has 
knowledge  of  only  the  overall  mean  period  demand  and  variance  in  the 
cycle.  Thus,  we  apply  the  power  approximation  policies  used  in  the 
stationary  model  to  the  nonstationary  models.  We  will  refer  to  these 
policies  as  stationary  policies. 

Table  2.19  contains  the  components  of  aggregate  total  co6t  per 
period  for  the  nonstationary  multi-item  systems  when  controlled  by 
stationary  policies.  Table  2.20  displays  the  absolute  and  percentage 
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Table  2.19 

Avera  .±  Costs  Per  Period  for  Multi-Item  Negative  Binomial  Systems: 
Controlled  With  Stationary  Policies 


COST  COMPt  ENT 

STATIONARY 

MODEL 

MODEL 

1 

MODEL 

II 

MODEL 

III 

MODEL 

IV 

INVENTORY 

440  (57.8) 

451  (45.7) 

447  (48.9) 

492  (22.1) 

466  (35.6) 

BACKLOG 

109  (14.3) 

326  (33.1) 

258  (28.2) 

1530  (68.8) 

636  (48.6) 

REPLENISHMENT 

212  (27.8) 

209  (21.2) 

210  (23.0) 

201  (9.0) 

206  (15.8) 

TOTAL 

761  (100.0) 

986  (100.0) 

915  (100.0) 

2222  (100.0) 

1308  (100.0) 

Note:  Figures  in  parentheses  are  percent  of  total. 


Table  2.20 

Average  Costs  Per  Period  for  Multi-Item  Negative  Binomial  Systems: 
Comparison  of  Stationary  Control  with  Optimal 
Control  Given  Full  Information 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

II 

MODEL 

IV 

INVENTORY 

1 (.1) 

18  (4.1) 

10  (2.2) 

63  (14.6) 

32  (7.4) 

BACKLOG 

0 (-.2) 

221  (208.7) 

151  (142.0) 

1428  (1406.1) 

532  (510.9) 

REPLENISHMENT 

1 (.3) 

1 (.4) 

2 (1.2) 

13  (6.9) 

5 (2.4) 

TOTAL 

1 (.1) 

239  (32.0) 

163  (21.7) 

1504  (209.4) 

i 

569  (77.0) 

Note:  Table  2.20  shows  the  absolute  increases  in  the  cost  components  of  Table 
2.19  over  those  of  Table  2.2  with  percentage  increases  shown  in 
parentheses . 
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Incremental  changes  in  costs  when  stationary  policies  are  compared 
to  optimal  control  given  full  information. 

Table  2.20  clearly  indicates  that  stationary  policies  used  in 
nonstat lonary  environments  incur  large  percentage  increases  in  expected 
total  cost.  Model  1 has  the  smallest  percentage  increase  (32%)  in 
expected  total  cost.  Models  I and  II,  however,  have  percentage  in- 
creases in  total  cost  that  are  much  smaller  than  those  of  models  III 
and  IV.  This  suggests  that  larger  variations  in  demand  means  contri- 
butes to  the  poor  performance  of  the  stationary  policies.  The  increases 
in  backlog  cost  component  increases,  varying  from  142%  to  1406%,  are 
the  primary  source  of  excessive  costs  in  all  the  nonstationary  systems. 
With  the  possible  exception  of  model  III,  replenishment  and  inventory 
costs  do  not  differ  greatly  from  those  obtained  under  optimal  control. 

The  distribution  of  expected  average  costs  for  the  various  para- 
meter settings  is  shown  in  Table  2.21.  The  large  penalty  costs  of  the 
nonstationary  systems  result  in  different  apportions  for  the  components 
of  total  cost  when  contrasted  with  the  stationary  system.  The  total 
cost  distribution  of  model  I and  model  II,  however,  is  similar  to  that 
of  the  stationary  system,  while  that  of  model  III  and  model  IV  bears 
little  resemblance  to  the  stationary  system.  (Recall  in  Section  2.4 
that  under  approximately  optimal  control  the  distribution  of  total  cost 
for  all  systems  was  similar  to  that  of  the  stationary  model.)  This 
again  suggests  that  greater  variation  in  demand  means  results  in  further 
discrepancies  in  system  characteristics  from  the  stationary  model  under 
stationary  control. 

Table  2.22  records  the  percentage  increases  over  optimal  control  for 
each  component  of  cost  for  the  various  parameter  settings.  For  all  the 
nonstationary  systems  the  backlog  component  is  the  primary  source  of  ex- 
cessive cost.  This  is  especially  apparent  for  those  items  with  unit 
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penalty  cost  of  99.  Percentage  Increases  In  Inventory  and  replenish- 
ment costs  are  small  in  comparison  to  the  backlog  component  for  all 
parameter  settings  in  the  nonstationary  systems. 

The  values  for  selected  operating  characteristics  are  displayed 
in  Table  2.23.  Table  2.24  gives  the  percentage  change  in  the  value  of 
these  operating  characteristics  from  the  corresponding  ones  in  Table  2.4. 
Backlog  frequencies  are  extremely  high  for  all  nonstationary  systems. 

For  those  items  with  unit  costs  of  99,  the  percentage  increases  range 
from  156.2%  in  model  II  to  1033.6%  in  model  III. 

In  summary,  stationary  policies  are  unacceptable  approximations  in 
a nonstationary  environment.  The  desirable  simplicity  of  stationary 
policies  does  not  adequately  compensate  for  the  excess  values  of  backlog- 
related  operating  characteristics. 

We  now  examine  the  situation  where  the  decision  maker  has  knowledge 
of  the  basic  movement  of  demand  means  through  the  cycle,  but  incorrect 
information  about  the  absolute  values  of  the  demand  means.  Model  III 
is  selected  for  a detailed  examination  under  this  type  of  misspecif ication. 

We  assume  that  we  correctly  identify  the  true  values  of  the  demand 
means  for  the  first  seven  periods  of  the  cycle.  During  the  last  five 
periods  of  the  cycle,  however,  we  assume  that  the  mean  demand  gradually 
Increases  to  only  twice  the  level  of  the  first  seven  periods,  reaching 
a maximum  value  in  period  10,  and  then  gradually  decreases.  In  other 
words,  we  expect  mean  demand  to  peak  at  a much  lower  mean  demand  level 
than  is  actually  the  case  for  model  III.  The  variance- to-mean  ratio  is 
still  assumed  constant  at  3,  so  that  a misspecif ication  in  the  mean  de- 
mand results  in  a misspecif ication  in  the  demand  variance.  We  now  inves- 
tigate the  consequences  of  such  a misspecif ication  on  the  system's  opera- 
ting characteristics  for  model  III.  The  power  approximation  policy 
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Percentage  Excess  of  Operating  Characteristics  of  a Multi-item  System  Under  Stationary 
Control  Over  Those  for  System  Under  Optimal  Control  with  Full  Information 
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resulting  from  such  a misspecif ication  will  be  referred  to  as  a 
misspecif ied  policy. 

Table  2.25  contains  the  aggregate  values  of  total  cost  per  period 
and  its  components  for  model  III  when  controlled  by  the  misspecif ied 
policies.  Backlog  cost  is  clearly  the  dominant  cost  component,  account- 
ing for  89.2%  of  the  total  cost.  Table  2.26  provides  comparisons  for 
the  costs  of  misspecif ied  policies  with  those  of  optimal  control.  Ap- 
parently the  smaller  misspecified  peaks  in  demand  result  in  insufficient 
levels  of  inventory.  Holding  costs  are  54.9%  below  their  optimal  levels. 

As  a result,  backlog  costs  are  over  3400%  in  excess  of  optimal  control. 

This  huge  backlog  cost  is  chiefly  responsibile  for  a percentage  increase 
in  total  cost  of  462.2%  for  the  multi-item  system. 

More  detailed  breakdowns  of  Tables  2.25  and  2.26  are  provided  in 
Tables  2.27  and  2.28.  Although  large  percentages  of  excess  backlog 
cost  are  found  under  every  parameter  setting,  large  unit  penalty  costs, 
leadtimes,  or  means  result  in  the  greatest  increment  in  backlog  costs. 

In  contrast,  the  percentage  decreases  of  inventory  costs  are  much  more 
uniform  over  all  of  the  parameter  settings. 

The  values  of  operating  characteristics  are  contained  in  Table 
2.29  with  the  percentage  excess  over  optimal  control  recorded  in  Table  2.30. 
The  higher  values  of  backlog  cost  noted  under  large  unit  penalty  costs, 
leadtimes,  or  means  are  also  reflected  in  the  corresponding  values  of 
backlog  frequency  and  weighted  proportion  of  demand  backlogged  in 
Table  2.29. 

Hence,  we  conclude  that  the  savings  in  inventory  costs  that  result 
from  the  misspecified  demand  structure  are  not  sufficient  to  compensate 
for  the  increase  in  backlog  costs  and  backlog  frequencies.  Precise 
knowledge  of  the  movement  of  demand  means  and  variances  throughout  the 
cycle  is  crucial  in  obtaining  approximately  optimal  policies. 
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Table  2.25 
Ml 

[ Average  Costs  Per  Period  for  a 

| Multi-Item  Negative  Binomial  System: 

Controlled  With  Misspecified  Policies 


COST  COMPONENT 

MODEL 

III 

INVENTORY 

193 

(4.8) 

BACKLOG 

3606 

(89.2) 

REPLENISHMENT 

243 

(6.0) 

TOTAL 

4042 

(100.0) 

Note:  Figures  in  parentheses  are  percent  of  total. 


Table  c.26 


» Average  Costs  Per  Period  for  a Multi-Item 

Negative  Binomial  System:  Comparison  of 
Misspecified  Control  with  Optimal 
Control  Given  Full  Information 


COST  COMPONENT 

MODEL 

III 

INVENTORY 

-236 

(-54.9) 

BACKLOG 

3505 

(3450.8) 

REPLENISHMENT 

5555 

(29.3) 

TOTAL 

3324 

(462.9) 

Note:  Table  2.26  shows  the  absolute  Increases  in  the 
cost  components  of  Table  2.25  over  those  of 
Table  2.2  with  percentage  increases  shown  in 
parentheses . 


Percentage  Excess  of  Costs  of  a Multi-Item  System  Under  Misspecified  Control  Over  the  Costs  of 
Controlling  the  Same  System  Optimally  with  Full  Information 
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Operating  Characteristics  of  a Multi-Item  System  Controlled  with 
Misspecified  Information  About  Demand 
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3.  SYSTEM  CONTROL  WITH  STATISTICAL  INFORMATION  ABOUT  DEMAND 

We  suppose  that  a manager  responsible  for  the  control  of  an  inven- 
tory system  as  defined  in  this  section  knows  the  cycle  length  and  the 
pattern  of  demand  means  through  the  cycle,  that  is,  each  mean  relative 
to  the  smallest  demand  mean.  It  was  shown  in  Section  1.2.3.  that  esti- 
mates of  the  mean  and  variance  of  the  monthly  demand  distributions  can 
be  obtained  from  this  information  and  an  historical  sample  of  realized 
demand.  Such  estimates  are  used  below  instead  of  the  values  for  the 
true  means  and  variances  in  control  policies  that  are  approximately 
optimal  in  the  case  of  perfect  information.  Details  of  the  control 
policies  were  given  previously  in  Section  1.1.2. 

The  operating  cnaracteristics  describing  the  performance  of  inven- 
tory systems  under  control  using  only  statistical  information  about 
demand  have  not  been  computed  analytically.  The  results  in  this  section 
are  statistical  estimates  of  the  performance  of  some  selected  inventory 
systems  found  by  computer  simulation.  During  each  simulation  run  the 
control  policy  parameters  are  revised  at  regular  intervals  to  allow  for 
potential  changes  in  the  underlying  demand  structure.  Such  changes  do 
not,  in  fact,  occur  during  the  simulation  experiments,  but  in  a real- 
life  system,  periodic  revision  is  essential  to  detect  and  adapt  to  a 
possibly  changing  demand  structure. 

Section  3.1  discusses  the  results  from  the  simulation  of  multi- 
item systems  under  statistical  control.  Revision  intervals  and  history 
lengths  are  set  to  24  periods,  corresponding  to  two  years  of  previous 
data.  The  sensitivity  of  the  operating  characteristics  to  values  of 
the  input  parameters  is  Investigated. 

Under  statistical  control,  the  expected  cost  per  period  of  the 
stationary  system  is  estimated  to  be  7.1%  higher  than  when  the  system 
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is  controlled  optimally  with  full  information.  For  the  nonstationary 
systems,  the  estimated  percentage  increase  in  expected  total  cost  per 
period  depends  on  the  nature  of  the  nonstationarity , varying  from  7.9% 
to  14.8%. 

3 . 1 Performance  under  Control  with  Statistical  Information  about  Demand 

The  input  parameters  for  each  item  in  the  stationary  and  four 
nonstationary  16-item  systems  were  specified  in  Section  1.2.2.  Esti- 
mated average  total  cost  per  period  and  its  components  are  shown  in 
Table  3.1.  Under  statistical  control,  holding  and  backlog  costs  are 
significantly  higher  in  each  system  than  those  found  under  optimal 
control,  shown  in  Table  2.2.  With  the  exception  of  model  III,  however, 
replenishment  costs  are  almost  identical. 

Table  3.2  shows  estimates  of  absolute  and  percentage  incremental 
changes  in  expected  costs  per  period  when  statistical  control  is  com- 
pared with  optimal  control  given  full  information  (cf.  Table  2.2).  The 
stationary  model  has  the  smallest  percent  increase  in  expected  total 
cos t , 7.1%.  Model  I and  model  II  have  increases  only  slightly  above  that 
amount,  but  model  III  and  model  IV  have  percentage  Increases  of  14.8% 
and  10.6%,  respectively.  Apparently  larger  variations  in  demand  means 
produce  greater  percentage  increases  in  total  cost  for  model  III  and 
model  IV.  All  systems  studied  incurred  excess  inventory  and  backlog 
costs. 

When  the  items  in  each  system  are  divided  into  subsystems  accord- 
ing to  values  taken  by  input  parameters,  the  resulting  apportionment  of 
costs  between  the  subsystems  is  shown  in  Table  3.3.  The  distribution  of 
costs  are  similar  to  those  in  Table  2.3,  which  shows  the  same  systems 
controlled  with  full  information. 
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Tabic  3.1 

Average  Cotta  Per  Period  for  Multi-Item  Negative  Binomial  Syateas: 

Controlled  With  Statistical  Information  About  Demand 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

11 

MODEL 

III 

MODEL 

IV 

INVENTORY 

463 

(56.9) 

473 

(58.3) 

469 

(57.8) 

496 

(60.2) 

492 

(60.1) 

BACKLOG 

138 

(16.9) 

130 

(16.1) 

133 

(16.4) 

128 

(15.6) 

124 

(15.2) 

REPLENISHMENT 

213 

(26.2) 

208 

(25.6) 

209 

(25.8) 

199 

(24.2) 

202 

(24.7) 



TOTAL 

813 

(100. 0) 

810 

(100.0) 

811 

(100.0) 

824 

(100.0) 

818 

(100.0) 

Note'  Figures  In  parentheses  are  percent  of  total. 


Table  3.2 

Average  Costs  Per  Period  for  Multi-Item  Negative  Binomial  Systems: 
Comparison  of  Statistical  Control  With  Optimal 
Control  Clven  Full  Information 

(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

11 

MODEL 

III 

MODEL 

IV 

INVENTORY 

24 

(5.4) 

40 

(9.2) 

31 

(7.2) 

68 

(15.8) 

58 

(13.4) 

BACKLOG 

28 

(25.9) 

■ 24 

(23.2) 

26 

(25.0) 

27 

(26.5) 

20 

(19.2) 

REPLENISHMENT 

2 

(.9) 

0 

(-.2) 

1 

(.7) 

12 

(6.2) 

I 

(.3) 

TOTAL 

54 

(7.1) 

64 

(8.6) 

59 

(7.9) 

106 

(14.8) 

79 

(10.6) 

Note:  Table  3.2  shows  the  absolute  Increases  In  the  cost  components  of  Table  3.1 
over  those  of  Table  2.2  with  percentage  Increases  shown  In  parentheses. 


Apportionment  of  Costs  for  a Multi-Item  System  Controlled  with  Statistical  Information  About  Demand 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 
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Table  3.4  shows  the  estimated  percentage  Incremental  changes  in  the 
expected  average  total  cost  and  Its  components.  For  the  various  sub- 
systems, this  table  elaborates  on  the  changes  summarized  before  In  Table 
3.2. 

An  examination  of  the  percentage  Increases  indicates  that  for  all 
systems  the  excess  backlog  costs  noted  In  Table  3.2  are  due  mainly  to 
those  items  with  high  unit  penalty  costs.  This  is  also  true  to  a lesser 
extent  for  inventory  costs.  It  should  be  noted,  however,  that  the  values 
in  Table  3.4  are  only  estimates  of  the  increase  in  costs  of  controlling 
the  systems  when  the  control  is  statistical  instead  of  with  full  infor- 
mation, as  rules  for  optimal  control  in  the  identical  circumstances  of 
statistical  information  and  periodic  revision  have  not  been  derived 
[cf.  Hayes  (1969)]. 

I The  physical  operating  characteristics  of  the  multi-item  system  have 

their  estimated  values  tabulated  in  Table  3.5,  and  Table  3.6  sets  out  the 
percentage  increase  of  each  over  the  corresponding  values  found  in  Table 
2.4.  As  in  the  case  of  approximately  optimal  control,  model  III  and 
model  IV  have  higher  inventory  levels  and  lower  backlog  frequencies  than 
the  other  systems. 

The  direction  of  the  sensitivity  of  the  policy  variables  and  the 
operating  characteristics  to  changes  in  the  input  parameter  is  indicated 
in  Table  3.7.  This  table  shows  that  the  behavior  of  the  stationary  and 
nonstationary  systems  has  the  same  pattern  with  one  exception.  Whereas 
the  backlog  frequency  of  model  III  is  monotonically  increasing  with 
respect  to  mean  demand  and  monotonically  decreasing  with  respect  to 
replenishment  cost,  the  reverse  is  true  for  the  other  systems. 


Some  differences  are  apparent,  however,  when  behavior  is  compared  to 
the  system  under  optimal  control  (Table  2.5).  The  monotonic  patterns,  for 
example,  of  expected  replenishment  cost,  expected  replenishment  frequency, 


Table  3.5a 


WEIGHTED  PROPORTION  OF 

DEMAND  BACKLOGGED  .013  .160  .007  .013  .014  .012  .015  .017 

REPLENISHMENT  FREQUENCY  .283  .286  .281  -322  .244  .288  .278  .236 


Operating  Characteristics  of  a Multi-Item  System  Controlled  with  Statistical  Information  About  Demand 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


REPLENISHMENT  FREQUENCY  .286  .288  .284  .325  .246  .289  .282  .238  .334 


Table  3.5c 
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Operating  Characteristics  of  a hulti-Item  System  Controlled  with  Statistical  Information  About 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


Table  3.6a 


Percentage  Excess  of  Operating  Characteristics  of  a Mult.'-Item  System  Under  Statistical 
Control  Over  Those  for  System  Under  Optimal  Contro.  with  Full  Information. 


L 2 2 


REPLENISHMENT  FREQUENCY  | .4  I 6.4  -5.0 


Table  3.6c 


Percentage  Excess  of  Operating  Characteristics  of  a Multi-Item  System  Under  Statistical 
Control  Over  Those  for  System  Under  Optimal  Control  with  Full  Information. 


2 


REPLENISHMENT  FREQUENCY  .2  6.3  -5. A 
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Table  3.7a 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of  an 
Inventory  System  Controlled  With  Statistical  Information  About  Demand 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT)/C(IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 

STATIONARY  MODEL 


OPERATING  CHARACTERISTIC 

E(COST)  t 

E (PERIOD- END  INVENTORY)  f 

E( BACKLOG  COST)  + 

E( BACKLOG  FREQUENCY  + 

PROPORTION  DEMAND  BACKLOGGED  + 
E (REPLENISHMENT  COST)  1 

E (REPLENISHMENT  FREQUENCY  + 

POLICY  PARAMETERS 

D -v 

s t 

S f 


f 

+ 

+ 

+ 

+ 

t 


+ 

+ 


1 

+ 

f 

+ 

x 


+ 

+ 


MODEL  I 

OPERATING  CHARACTERISTIC 

E(COST)  t + 

E (PERIOD- END  INVENTORY  t t 

E (BACKLOG  COST)  + t 

E( BACKLOG  FREQUENCY)  + + 

PROPORTION  DEMAND  BACKLOGGED  + t 

E (REPLENISHMENT  COST)  + t 

E (REPLENISHMENT  FREQUENCY)  i * 


+ 

+ 

t 

f 

t 

+ 

i 


POLICY  PARAMETERS 
D 


<\» 


7 


t *■  t 


S 


t t 


t indicates 
+ indicates 
i Indicates 
7 indicates 


relation  Is  monotone  Increasing 
relation  Is  monotone  decreasing 
relation  Is  approximately  constant 
no  conclusion  can  be  drawn  from  the  data 


t 


f 

+ 

t 

+ 

+ 

t 


1 

+ 

+ 


+ 

t 

+ 

t 

t 


+ 

+ 


a. 
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Table  3.7b 

Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of  an 
Inventory  System  Controlled  With  Statistical  Information  About  Demand 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT)/C(IN)  C (FIX) /C (IN ) LEADTIME  MEAN  DEMAND 

STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 


E(COST)  4 4 
E (PERIOD-END  INVENTORY)  4 t 
E( BACKLOG  COST)  4 4 
E (BACKLOG  FREQUENCY)  4-  4 
PROPORTION  DEMAND  BACKLOGGED  4 4 
E (REPLENISHMENT  COST)  4 4 
E( REPLENISHMENT  FREQUENCY)  4 4 


A 

4 

4 

4 

4 

4 

4 


POLICY  PARAMETERS 

\ D -v  4 


s 4 4 

S 4 4 

MODEL  II 

OPERATING  CHARACTERISTIC 

E(COST)  4 4 

E (PERIOD- END  INVENTORY)  4 4 

E (BACKLOG  COST)  4 4 

E( BACKLOG  FREQUENCY)  4 4 

PROPORTION  DEMAND  BACKLOGGED  4 4 

E (REPLENISHMENT  COST)  4 4 

E (REPLENISHMENT  FREQUENCY)  4 4 

POLICY  PARAMETERS 

D -v  4 


4 

4 

4 


4 

4 

4 

4 

4 

4 

4 

? 


a 4 

S 4 

4 indicates  relation  is  monotone  increasing 
+ indicates  relation  is  monotone  decreasing 
'v  indicates  relation  is  approximately  constant 
? indicates  no  conclusion  can  be  drawn  from  the  data 


4 4 

4 4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 


4 

4 

4 

4 

4 

4 

t 


4 

4 

4 


Table  3.7c 
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Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of  an 
Inventory  System  Controlled  With  Statistical  Information  About  Demand 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT) /C (IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 


STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 

E(COST)  4 

E (PERIOD-END  INVENTORY)  4 

E( BACKLOG  COST)  4 

E( BACKLOG  FREQUENCY)  4 

PROPORTION  DEMAND  BACKLOGGED  4 

E (REPLENISHMENT  COST)  4 

E (REPLENISHMENT  FREQUENCY)  4 

POLICY  PARAMETERS 

D *v< 

s 4 

S 4 

MODEL  III 

OPERATING  CHARACTERISTIC 

E(COST)  4 

E(PERIOD-END  INVENTORY)  4 

E (BACKLOG  COST)  4 

E( BACKLOG  FREQUENCY)  4 

PROPORTION  DEMAND  BACKLOGGED  4 

E (REPLENISHMENT  COST)  4 

E (REPLENISHMENT  FREQUENCY)  4 

POLICY  PARAMETERS 

D ^ 

s 4 

S 4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


4 4 4 
4 4 4 
4 4 4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


4 

4 

4 


? 

4 

4 


* 

4 

4 


4 indicates  relation  is  monotone  increasing 
4 indicates  relation  is  monotone  decreasing 
^ indicates  relation  is  approximately  constant 
7 indicates  no  conclusion  can  be  drawn  from  the  data 
* indicates  an  explanation  is  given  in  the  text 
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Table  3. 7d 


Sensitivity  of  Input  Parameters  of  Operating  Characteristics  of  an 
Inventory  System  Controlled  With  Statistical  Information  About  Demand 
When  Demand  Has  a Negative  Binomial  Distribution 


INPUT  PARAMETERS 

C(OUT)/C(IN)  C(FIX)/C(IN)  LEADTIME  MEAN  DEMAND 

STATIONARY  MODEL 

OPERATING  CHARACTERISTIC 


E(COST)  + 
E (PERIOD- END  INVENTORY)  t 
E(BACKLOC  COST)  t 
E (BACKLOG  FREQUENCY)  + 
PROPORTION  DEMAND  BACKLOGGED  + 
E( REPLENISHMENT  COST)  + 


E (REPLENISHMENT  FREQUENCY)  + 

POLICY  PARAMETERS 


s t 

S t 

MODEL  IV 

OPERATING  CHARACTERISTIC 

E(COST)  t 

E (PERIOD-END  INVENTORY)  t 

E (BACKLOG  COST)  t 

E (BACKLOG  FREQUENCY)  + 

PROPORTION  DEMAND  BACKLOGGED  + 
E (REPLENISHMENT  COST)  + 


E (REPLENISHMENT  FREQUENCY)  4- 


POLICY  PARAMETERS 
D 
8 
s 


'X, 

t 

t 


4 

t 

+ 

+ 

t 

i 

4 


4 

4 

4 


4 

4 

4 

4 

4 

4 

4 


t 

4 

4 


t 

4 

4 

4 

4 

4 

4 

t 

4 

4 


4 

4 

4 

4 

4 

4 

i 


? 

4 

4 


t indicates 
4 indicates 
^ indicates 
? indicates 


relation  is  monotone  increasing 
relation  is  monotone  decreasing 
relation  is  approximately  constant 
no  conclusion  can  be  drawn  from  the  data 


t 

t 

4 

4 

4 

4 

4 


t 

4 

4 


4 

4 

4 

4 

4 

4 

4 


t 

4 

4 


* 
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and  expected  backlog  costs  with  respect  to  penalty  cost  are  reversed  from 
those  observed  for  the  same  operating  characteristics  under  optimal  control. 

Policy  parameters  s,  S,  and  D appear  to  follow  the  same  relationships 
as  those  found  under  approximately  optimal  control.  The  occasional  depart- 
ure from  increasing  monotonicity  of  D with  respect  to  mean  demand  that  was 
noted  in  model  III  under  approximately  optimal  control  is  also  apparent 
under  statistical  control  in  Appendix  D. 

We  also  consider  the  value  taken  by  the  ratio  1/(1  + tr),  where  it 
is  the  unit  backlog  penalty  cost,  as  an  approximation  to  the  frequency  of 
backlog  thac  can  be  observed  for  a typical  item  in  a multi-item  system. 

Table  3.8,  which  is  derived  from  Table  3.5,  shows  that  when  ir  - 4, 

1/(1  + it)  is  a good  backlog  frequency  approximation  for  all  systems 
except  model  III  and  model  IV,  where  1/(1  + it)  is  a loose  upper  bound; 

1/ (1  + it)  is  a good  approximation  of  backlog  frequency  for  all  systems 
when  it  * 99,  although  it  is  slightly  less  than  the  backlog  frequencies 
for  each  system. 

Table  3.8 


Backlog  Frequencies  for  Multi-Item  Systems 
With  Negative  Binomial  Demands 

C(OUT)/C(IN)  (-*) 


4 

99 

1/U+n) 

.2 

.01 

SYSTEM: 

Stationary  Model 

.195 

.013 

Model  I 

.190 

.012 

Model  II 

.192 

.012 

Model  III 

.179 

.012 

Model  IV 

.178 

.011 
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3 . 2 Statistical  Control;  Single-Item  Systems 

The  average  (s^.S^)  values,  for  1 - 1,...,12,  of  the  statistical 

policy  are  recorded  for  case  A and  case  B in  Tables  3.9  and  3.10, 

respectively.  The  values  for  s^  do  not  differ  greatly  from  these 

obtained  under  approximately  optimal  control.  The  statistical  values 

for  Sj  are  also  close  to  those  obtained  under  approximately  optimal 

— (3) 

control  with  full  information  with  one  exception:  when  Is  at  Its 

largest  value,  the  corresponding  S^  of  the  statistical  policy  is 

significantly  less  than  the  approximately  optimal  S^  . This  occurs 

in  period  10  in  the  cycle  for  case  A and  periods  8 and  9 for  case  B. 

— (3) 

The  statistical  estimates  for  during  these  periods  are  occasion- 

ally much  greater  than  the  actual  parameters.  As  mentioned  previously 
in  Section  2.4,  large  means  can  result  in  small  values  for  - s^f 
that  may  cause  the  power  approximation  to  produce  small  values  for  S^. 
This  explains  the  discrepancies  in  under  statistical  and  approxi- 

mately optimal  control.  This  difference,  however,  does  not  cause  the 
general  shape  of  the  cost  components,  displayed  in  Figure  3.1  for 
case  A and  Figure  3.2  for  case  B,  to  differ  from  those  obtained  under 
approximately  optimal  control,  shown  in  Figures  2.4  and  2.5. 


Table  3.9 

Statistical  Policy  for  Case  A 


Period  in 
C 


1 

58.4 

(57) 

[57] 

2 

58.4 

(57) 

[57] 

3 

58.4 

(57) 

[57] 

4 

58.4 

(57) 

[57] 

' 5 

58.4 

(57) 

[57] 

6 

58.4 

(57) 

[59] 

7 

58.4 

(57) 

[55] 

8 

73.9 

(72) 

[72] 

9 

89.1 

(87) 

[86] 

10 

104.0 

(102) 

[105] 

11 

89.1 

(87) 

[89] 

12 

73.9 

(72) 

[73] 

86.5 

(85) 

[82] 

86.5 

(85) 

[83] 

86.5 

(85) 

[82] 

86.5 

(85) 

[82] 

86.5 

(85) 

[84] 

86.5 

(85) 

[77] 

86.5 

(85) 

[67] 

105.8 

(104) 

[108] 

123.6 

(122) 

[130] 

131.2 

(140) 

[128] 

123.6 

(122) 

[112] 

105.8 

(104) 

[95] 

Note:  The  approximately  optimal  policy  Is  given  In  parentheses 
and  the  optimal  policy  is  contained  in  brackets. 


TGTflL  COST 


Components  for  Case  A Controlled  With  Statistical  Information  About  Demand 


Period  in 


Table  3.10 

Statistical  Policy  for  Case  B 


46.4 

(45) 

[46] 

35.4 

(35) 

[33] 

35.4 

(35) 

[33] 

46.4 

(45) 

[44] 

62.2 

(61) 

[60] 

76.6 

(76) 

[75] 

92.7 

(91) 

[90] 

102.6 

(101) 

[101] 

102.6 

(101) 

[102] 

92.7 

(91) 

[93] 

76.6 

(76) 

[78] 

62.2 

(61) 

[62] 

71.7 

(69) 

[61] 

56.5 

(56) 

[55] 

56.5 

(56) 

[56] 

71.1 

(69) 

[73] 

91.3 

(90) 

[93] 

110.4 

(109) 

[113] 

127.3 

(127) 

[128] 

132.0 

(139) 

[136] 

132.0 

(139) 

[131] 

127.3 

(127) 

[115] 

110.4 

(109) 

[95] 

91.3 

(90) 

[76] 

Note:  The  approximately  optimal  policy  is  given  in  parentheses 

and  the  optimal  policy  is  contained  in  brackets. 
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Components  for  Case  B Controlled  With  Statistical  Information  About  Demand 


4.  FORECASTING  THE  OPERATING  CHARACTERISTICS  OF  INVENTORY 
SYSTEMS  UNDER  STATISTICAL  CONTROL 


The  behavior  of  an  inventory  system  controlled  with  statistical  in- 
formation about  demand  was  investigated  in  Section  3 using  simulation 
techniques.  In  the  absence  of  analytic  expressions  for  the  operating 
characteristics  as  functions  of  the  input  distribution  functions  and 
other  system  parameters,  sample  distributions  for  the  operating  charac- 
teristics were  gathered  from  the  simulations,  particular  attention  being 
paid  to  the  estimation  of  expected  values. 

The  manager  of  an  inventory  system  operating  with  incomplete  infor- 
mation about  demand  will  need  to  forecast  the  system's  behavior,  especially 
to  justify  the  installation  of  scientific  control.  Forecasting  also  may 
be  a routine  requirement  consequent  upon  the  periodic  revision  of  the  para- 
meters of  the  replenishment  policies  or  even  as  part  of  a regular  budgeting 
procedure . 

Forecasting  methods  have  been  considered  in  MacCormick  [(1974), 

Section  3.2].  In  Section  4 information  is  gathered  about  forecasts  using 
the  same  techniques  as  MacCormick  which  involve  double  use  of  the  demand 
history,  for  control  and  forecasting.  Sample  distributions  of  the  fore- 
casts are  compared  with  those  of  the  operating  characteristics,  which  were 
obtained  in  Section  3. 

The  major  issues  to  be  resolved  here  are  the  extent  of  bias  and  the 
level  of  dispersion  of  the  forecasts,  together  with  the  variation  of 
these  with  system  settings,  including  the  demand  process,  unit  costs,  and 
the  replenishment  leadtime.  Little  consideration  is  given  to  an  Important 
sampling  problem,  namely,  that  since  it  is  usually  prohibitively  expensive 
to  gather  a history  of  demands  for  every  item  in  a large  multi-item  system. 
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a sample  of  the  items  must  be  selected  to  make  a forecast.  The  accuracy 
of  such  a forecast  Is  likely  to  be  critically  dependent  on  the  number  of 
items  in  the  sample  and  its  composition.  Forecasts  made  for  the  multi- 
item system  in  this  simulation  experiment  have  used  every  item  in  the 
system. 

4 . 1 Properties  of  Forecasts  for  Operating  Characteristics  of  Multi-Item 
Systems 

In  the  simulation  experiments,  point  forecasts  of  the  expected  values 
of  the  system  operating  characteristics  are  made  each  time  the  control 
parameters  (s,S)  are  revised.  The  demand  history  used  to  make  a forecast 
has  been  selected  to  contain  demands  for  the  same  number  of  periods  as  the 
revision  interval.  The  actual  values  for  the  forecasts  are  obtained  by 
retrospective  simulation,  which  has  been  explained  previously  in  Section 
1.2.6. 

Sample  average  values  of  the  forecasts  of  average  costs  per  period 
for  a multi-item  system  are  shown  in  Table  4.1,  the  sample  standard  de- 
viations being  shown  in  parentheses.  The  sample  statistics  for  forecasts 
shown  in  Table  4.1  can  be  compared  with  those  in  Table  3.1  for  the  average 
cost  per  period.  The  comparison  for  expected  average  cost  per  period  is 
explicit  in  Table  4.2,  which  contains  the  percentage  difference  between 
the  forecasts  and  the  estimated  values. 

On  the  average,  every  component  of  expected  cost  is  underestimated 
by  its  forecast,  behavior  predicted  by  the  analysis  of  simplified  models 
in  Section  2 of  MacCormick  (1974).  For  expected  total  cost  per  period, 
the  underestimation  is  as  high  as  10.8%  in  the  case  of  model  111.  For  all 
systems  the  bias  is  greatest  for  the  backlog  component. 

To  check  the  sensitivity  of  the  bias  of  the  forecasts  to  other  sys- 
tem parameters,  Table  4.3  was  constructed  showing  percentage  values  for 
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Table  4.1 


Avarag*  Forecasted  Values  for  txpactei  Coata  Par  Prriod  for  Mulrl-Ii.m  Negative  »lno«lal  Sy.t.-au. 
Controlled  With  )tatlatlcal  Information  About  Oaoumd 
(Revision  Intarval  2a  Periods.  Revision  lllatory  Length  2a  Parloda) 


COST  COMPONENT 

STATIONARY 

MODEL 

MODEL 

I 

MODEL 

II 

MODEL 

III 

) 

MkiaiEl. 

iv 

1 

INVENTORY 

4 SO  (20)  [60. R] 

454  (21)  [61.6] 

452  (21)  [61.3] 

475  (28)  [64.6] 

I 470  (24)  [b  1.  b] 

BACKLOG 

78  (28)  [10.6] 

75  (26)  [10. l] 

76  (30)  [10.3] 

62  (27)  [8.4] 

a;  uo 

REPLENISHMENT 

212  (5)  [28.7] 

209  (6)  [28.3] 

209  (5)  [28.4] 

199  (6)  [27.0] 

202  (6)  [27.3] 

TOTAL 

740  (36)  [100. 0] 

738  (34)  [100. 0] 

738  (38)  [100.0] 

735  (38)  [100.0] 

739  (35)  [100.0] 

Notn:  Estimates  of  the  standard  deviations  for  the  forecasts  are  shown  In  parentheses, 
and  each  coat  component's  percent  of  the  total  la  shown  in  brackets. 


Table  4.2 

Average  Coats  Per  Period  for  Multi-Item  Negative  Binomial  Systems 
Controlled  With  Statistical  Information  About  Demand: 
Comparison  of  Forecasts  With  Corresponding  Estimated  Expected  Values 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


COST  COMPONENT 

STATIONARY  1 
MODEL 

MODEL 

I 

MODEL 

II 

MODEL 

III 

MODEL 

IV 

INVENTORY 

2.9 

4.0 

3.6 

4.3 

4.3 

BACKLOG 

43.1 

42.7 

42.7 

52.1 

46.1 

REPLENISHMENT 

.4 

-.6 

-0.0 

.2 

.1 

TOTAL 

9.0 

9.0 

9.0 

10.  B 

9.6 

Note:  Table  4.2  shows  the  percentage  value  of  the  underestimate  by  the 
fort  cast  of  the  realized  value. 
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Forecasting  for  a Multi-Item  System  Under  Statistical  Control: 

Estimated  Percentage  Difference  Between  Forecasts  and  Expected  Values  of  Operating  Characteristics 
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the  biases.  In  every  system  examined,  the  bias  in  the  forecast  of  expected 
average  total  cost  tends  to  increase  with  an  increase  in  the  unit  backlog 
cost  or  in  the  mean  demand.  A very  large  increase  in  bias  is  observed  in 
backlog  cost  when  unit  penalty  costs  increases  from  4 to  99.  Table  4.4 
demonstrates  biases  for  forecasts  of  operating  characteristics  other  than 
costs,  expressed  as  percentages  of  the  estimates  of  the  actual  expected 
values  of  the  operating  characteristics.  In  particular,  the  forecasts  of 
the  backlog  frequency  do  not  have  a percentage  bias  as  large  as  those  for 
the  weighted  proportion  of  demand  backlogged.  For  both  characteristics, 
however,  the  percent  bias  is  much  greater  for  those  items  with  unit  backlog 
costs  of  99  than  those  with  unit  backlog  costs  of  4.  More  detailed  versions 
of  Tables  4. '3  and  4.4  can  be  found  in  Appendices  M,  Q,  U,  and  Y. 

It  is  important  to  note  that  each  of  the  200  sample  values  for  an 
operating  characteristic  that  are  collected  during  the  simulation  of  a 
single  item  is  based  on  24  demand  periods  of  the  simulation.  Since  24  Is 
a multiple  of  the  demand  cycle  length,  which  is  12,  the  first  demand  en- 
countered in  collecting  each  of  the  200  sample  forecast  values  is  always 
sampled  from  the  same  demand  distribution  in  the  demand  cycle.  Demand 
distribution  12  - 2X  + 1,  where  X is  the  leadtime,  is  arbitrarily  se- 
lected by  the  simulation  as  the  initial  distribution.  Thus,  the  question 
arises  as  to  whether  initializing  with  a different  demand  distribution 
could  result  in  changes  in  the  forecasts,  or  whether  24  demand  periods  is 
sufficient  to  diminish  the  effects  that  the  specific  initial  distribu- 
tion may  have  on  the  forecasted  operating  characteristic. 

This  issue  is  investigated  by  an  additional  simulation  of  model  IV, 
under  statistical  control,  using  a different  initializing  distribution. 

For  each  item  in  the  system  a new  initial  distribution  is  chosen  such  that, 
of  all  the  distributions  within  a cycle,  the  new  distribution’s  mean 
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differs  the  most,  in  absolute  value,  from  the  previously  selected  initial 
distribution.  Hence,  the  demand  distributions  in  periods  3 and  10  are 
used  as  initial  distributions  for  those  cases  with  leadtimes  of  2 and  4, 
respectively.  The  resulting  forecast  values,  shown  in  Table  4.5,  are 
virtually  the  same  as  the  original  simulation  of  model  IV.  We,  thus,  con- 
clude that  24  demand  periods  is  sufficient  to  deter  any  potential  effect 
of  an  initial  distribution,  and  that  this  issue  does  not  warrant  any 
further  investigation. 

Table  4.5 

Average  Forecasted  Values  for  Expected  Costs  Per  Period  for 
Multi-Item  Negative  Binomial  Systems 
Controlled  With  Statistical  Information  About  Demand 
(Revision  Interval  24  Periods,  Revision  History  Length  24  Periods) 


COST  COMPONENT 

MODEL  IV 

(From  Table  4.1) 

MODEL  IV 

(Using  a different 
init ializing 
distribution) 

INVENTORY 

470 

(24) 

[63.6] 

472 

(24) 

[64.0] 

BACKLOG 

67 

(26) 

[9.0] 

65 

(26) 

[8.7] 

REPLENISHMENT 

202 

(6) 

[27.3] 

202 

(6) 

[27.3] 

TOTAL 

735 

(35) 

[100.0] 

741 

(36) 

[100.0] 

Note:  Estimates  of  the  standard  deviations  for  the  forecasts 

are  shown  in  parentheses,  and  each  cost  component's  percent 
of  the  total  is  shown  in  brackets. 


4 . 2 Properties  of  Forecasts  for  Operating  Characteristics  of  Single-Item 
Systems 

The  average  forecasted  values  of  cost  components  are  displayed  in 
Figures  4.1  and  4.2,  corresponding  to  cases  A and  B.  The  underestimai ion 
of  backlog  costs,  observed  In  Table  4.1,  for  the  multi-item  systems  be- 
comes apparent  in  cases  A and  B when  Figures  4.1  and  4.2  are  compared  with 
Figures  3.1  and  3.2.  Both  single-item  systems  considerably  underestimate 
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their  backlog  costs  Incurred  under  statistical  control. 

Case  A's  increased  backlog  cost  in  period  9 under  statistical  control 
is  easily  seen  in  Figure  3.1.  This  increase  is  predicted  to  some  extent 
in  the  backlog  component  of  Figure  4.1,  but  the  forecasts  also  predict  non- 
existent increased  backlog  costs  in  periods  2 and  3.  The  backlog  cost 
forecasts  for  case  A underestimate  actual  backlog  costs  and  do  not  accurately 
reflect  the  behavior  of  the  backlog  component  in  the  demand  cycle.  The 
general  trends,  however,  of  holding  and  replenishment  costs  under  statisti- 
cal control  are  accurately  reflected  in  the  forecasts. 

Case  B also  fails  to  accurately  predict  statistical  backlog  costs. 

The  sinusoidal  pattern  of  backlog  costs  under  statistical  control,  observed 
in  Figure  3.2,  is  not  detectable  in  the  forecasts.  The  forecasted  backlog 
values  also  underestimate,  in  general,  the  statistical  values.  As  In  case 
A,  however,  case  B does  accurately  forecast  the  general  sinusoidal  shape 
of  statistical  holding  and  replenishment  costs. 

Table  A7  of  Appendices  A to  E sets  out,  for  each  item  in  the  five 
systems,  the  bias  estimated  for  the  forecast  of  total  cost  per  period. 
Comparison  of  these  biases  with  the  standard  deviations  estimated  for 
the  forecasts  of  total  costs,  set  out  in  Table  A8,  leads  to  the  conclu- 
sion that  except  for  certain  cases  with  unit  backlog  cost  of  99,  the  bias 
is  small  relative  to  the  standard  deviation  of  the  forecast,  or,  alterna- 
tively, relative  to  the  mean  square  error  taken  about  the  expected  average 
total  cost  per  period.  For  a single  item,  then,  there  is  little  point  to 
making  a correction  to  each  forecast  of  the  long  run  value  of  expected 
average  total  cost  per  period.  The  dispersion  of  the  forecasts  over- 
whelms any  bias.  It  already  has  been  observed  that  correction  for  bias 
is  worthwhile  for  a multi-item  system  containing  a sufficiently  large 
number  of  items  with  independent  demand  processes.  According  to  the  law 
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of  large  numbers,  asymptotically  the  dispersion  becomes  small  relative  to 
the  bias  in  the  forecast  of  an  aggregate  operating  characteristic. 

Table  A7  also  contains  information  useful  to  evaluate  the  forecasts 
for  the  short  run.  The  table  has  been  constructed  by  analysis  of  the  time 
series  of  differences  between  the  forecasts  and  the  values  subsequently 
realized  for  the  operating  characteristics.  A convention  has  been  adopted 
that  the  bias  is  positive  if  the  realized  value  of  an  operating  character- 
istic exceeds  on  the  average  the  forecast  value.  To  test  whether  bias 
observed  in  any  operating  characteristic  is  significantly  positive  or 
negative,  a two-tail  large  sample  test  has  been  applied  at  the  .05  level 
of  significance,  assuming  that  the  sample  mean  differences  for  200  obser- 
vations are  normally  distributed.  The  standard  deviation  used  to  construct 
the  unit  normal  test  statistic  is  that  estimated  for  the  sequence  of  200 
differences.  Results  of  the  tests  are  shown  in  columns  (6)  to  (11)  of 
Table  A7.  Note  that  the  average  total  cost  per  period  has  bias  signifi- 
cantly positive  for  all  the  items  in  each  of  the  five  systems.  A sig- 
nificantly positive  bias  also  is  present  for  average  period-end  inventory 
and  backlog  quantity  in  most  cases.  Stockout  frequency  is  significantly 
positive  when  the  unit  penalty  cost  is  large. 

In  summary,  for  single-item  systems  where  the  forecasts  are  to  be 
used  as  predictions  for  long  run  expected  values,  the  forecasts  are  esti- 
mated to  have  standard  deviations  large  relative  to  the  expected  values. 

They  are,  therefore,  not  useful  to  predict  well  the  long  run  performance 
of  a single-item  system.  For  a multi-item  system,  bias  in  the  forecasts 
has  been  observed  to  be  quite  general.  Both  this  bias  and  the  dispersion 
of  the  forecasts  are  sensitive  to  the  characteristics  of  the  items  in  the 


system,  however,  and  therefore  the  magnitudes  depend  on  the  overall  compo- 
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sition  of  the  multi-item  system  under  consideration.  For  this  reason, 
when  forecasting  is  to  be  done  using  only  a sample  of  items  in  the  system, 
care  must  be  taken  to  select  a sample  representative  of  the  overall 
structure  of  the  multi-item  system. 
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5.  TOPICS  FOR  FURTHER  RESEARCH 

In  this  section  we  summarize  the  main  results  of  this  report  and 
suggest  topics  for  further  investigation. 

In  Section  2.3  we  examined,  in  detail,  the  optimal  policies  and 
corresponding  operating  characteristics  under  optimal  control  for  two 
typical  single-item  cases.  We  observed,  in  particular,  that  the  movement 
of  optimal  policies  through  the  demand  cycle  was  similar  to  that  of  the 
means  and  variances  of  the  A + 1 - fold  demand  distributions,  + ^ » 

for  i = 1,...,  12.  This  suggested  that  an  approximately  optimal  (s,S) 
inventory  policy  can  be  computed  in  the  beginning  of  the  i*"^  period  within 
the  cycle  as  a function  of  the  cost  parameters  and  the  demand  means  and 
variances  for  only  periods  i to  A + 1 + i. 

The  12  convoluted  means  and  variances,  + and  0^(A+1)  ^ for 

i = 1,...,  12,  were  used  in  the  Ehrhardt  power  approximation  to  compute 
12,  possibly  distinct,  (s,S)  policies,  (see  Section  1.1.2  for  details). 
These  policies  resulted  in  operating  characteristics  that  were  very  close 
to  that  of  under  optimal  control. 

When  the  12  individual  demand  means  and  variances  were  not  known, 
they  were  estimated  using  previous  demand  data.  The  technique  adopted 
under  statistical  control  assumed  knowledge  of  the  relative  magnitudes 
of  demand  means  within  the  cycle  by  representing  mean  demand  in  period  1 
as  a x^,  wiiere  u is  the  smallest  demand  mean  encountered  in  the  cycle 
and  the  x^  variables  are  assumed  known.  a is  estimated  by  regression 
techniques,  (see  Section  1.2.3  for  details).  The  results  of  Section  1 
indicated  that  this  statistical  approach  resulted  in  relatively  small 


increases  of  total  costs  over  that  of  optimal  control. 


We  realize  that  in  this  study  the  amount  of  variation  of  mean  de- 
mand for  the  various  nonatationary  systems  has  been  arbitrarily  selected. 
Models  I and  II  vary  mean  demand  from  a base  level  a to  2a,  while 
models  III  and  IV  vary  mean  demand  from  a to  5u.  Sections  2 and  3 
reveal  a slight  deterioration  in  the  performance  of  the  power  approxi- 
nation  for  models  III  and  IV.  Hence,  it  may  be  interesting  to  increase 
the  mean  demand  levels  for  these  models  to  much  higher  levels  than  5a, 
to  see  if  further  degradation  occurs  In  the  event  of  further  degrada- 
tion, the  individual  variances  (j^C'+l)  ^ for  ^ m 1,...,12,  might  be 

slightly  altered  to  account  for  the  amount  of  variation  of  mean  demand 
within  a cycle.  Further  research  with  nonstationary  models  that  have 
small  variations  in  mean  demand  might  also  be  considered  to  see  if  the 
stationary  policies  In  Section  2.6  improve  their  performance. 

The  study  has  required  that  all  demands  are  stocastically  independent 
of  each  other.  Hence,  knowledge  of  simply  the  demand  in  one  period  pro- 
vides p r>  additional  information  about  future  demands.  Correlated  demand 
sequences  can  provide  an  interesting  area  for  further  investigation. 

The  analysis  of  inventory  policies  when  demands  follow  a Markov  process 
remains  an  open  area  for  research. 

Under  statistical  control  in  Section  3 we  assumed  knowledge  of  the 
general  shape  of  the  demand  means  within  a cycle  via  the  x^'s,  (see 
Section  1.2.3  for  details).  Complete  knowledge  of  the  x^'s  indirectly 
gives  us  the  cycle  length.  Relaxing  these  assumptions  would  require 
that  new  estimation  schemes  for  the  x^'s  and  the  cycle  length,  based 
on  the  observed  demands,  be  developed. 
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Appendix  A 


Single-Item  Data 


Stationary  Model 

Summary  of  Data  for  16  Items  with  Negative  Binomial  Demand 
Distributions  (Variance/Mean  = 3)  Controlled  with: 

Optimal  Policies  (DP) 

Power  Approximation  (PA) 

Statistical  Power  Approximation 


Table  A1 

Average  Total  Cost 

A 1 

A3 

Period-End  inventory 

A 2 

A4 

Period-End  Backlog  as  Proportion  of  Mean 

Demand 

A3 

A5 

Frequency  of  Periods  with  Backlog 

A4 

A6 

Replenishment  Frequency 

A5 

A7 

Estimated  Bias  of  Forecast  of  Total  Cost 

A6 

A8 

Est limited  Standard  Deviation  of  Forecast 
Total  Cost 

of 

A7 

A9 

Values  for  ( s ^ , S ^ ) 

AH  to  A1 3 

A10 

Standard  Deviations  of  (s^,S^)  Values 

A14  to  A1 

Note:  For  corresponding  data  in  MacCormick  (1974),  see  tables 
of  the  same  number  in  Appendices  A of  those  reports. 


Appendix  B 
Single-item  Data 

Model  I 

Summary  of  Data  for  16  Items  with  Negative  Binomial  Demand 
Distributions  (Variance/Mean  = 3)  Controlled  with: 

Optimal  Policies  (DP) 

Power  Approximation  (PA) 

Statistical  Power  Approximation 


Table  A1 

Average  Total  Cost 

nil i£ 

B) 

A3 

Period-End  Inventory 

B2 

A4 

Period-End  Backlog  as  Proportion  of  Mean 

Demand 

BJ 

A5 

Frequency  of  Periods  with  Backlog 

B4 

A6 

Replenishment  Frequency 

B6 

A7 

Estimated  Bias  of  Forecast  of  Total  Cost 

Bb 

A8 

Estimated  Standard  Deviation  of  Forecast 
Total  Cost 

of 

B7 

A9 

Values  for 

B8 

A10 

Standard  Deviations  of  (s  ,S^)  Values 

B14 

Note:  For  corresponding  data  in  MacCormlck  (1974),  see  tables 
of  the  same  number  in  Appendices  A of  those  reports. 


Appendix  C 


I 


Single- Item  Data 


Model  II 

Summary  of  Data  for  16  Items  with  Negative  Binomial  Demand 
Distributions  (Variance/Mean  = 3)  Controlled  with: 

Optimal  Policies  (DP) 

Power  Approximation  (PA) 

Statistical  Power  Approximation 


ri'Ui: 

Table  A1 

Average  Total  Cost 

Cl 

A3 

Period-End  Inventory 

C2 

A4 

Period-End  Backlog  as  Proportion  of  Mean 

Demand 

C3 

A5 

Frequency  of  Periods  with  Backlog 

C4 

A6 

Replenishment  Frequency 

C3 

A7 

Estimated  Bias  of  Forecast  of  Total  Cost 

Cb 

A8 

Estimated  Standard  Deviation  of  Forecast 
Total  Cost 

of 

07 

A9 

Values  for  (s^.S^) 

C8  to  Cl J 

A10 

Standard  Deviations  of  (s.,S^)  Values 

Cl  4 to  Cl 

Note:  F >r  corresponding  data  in  MacCormick  (1974),  see  tables 

of  the  same  number  in  Appendices  A of  those  reports. 
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Appendix  D 


Single-item  Data 


Model  III 

Summary  of  Data  for  16  Items  with  Negative  Binomial  Demand 
Distributions  (Variance/Mean  = 3)  Controlled  with: 

Optimal  Policies  (DP) 

Power  Approximation  (PA) 

Statistical  Power  Approximation 


iil'il1: 

Table  A1 

Average  Total  Cost 

D1 

A3 

Period-End  Inventory 

D2 

A4 

Period-End  Backlog  as  Proportion  of  Mean 

Demand 

03 

A5 

Frequency  of  Periods  with  Backlog 

04 

A6 

Replenishment  Frequency 

D5 

A7 

Estimated  Bias  of  Forecast  of  Total  Cost 

06 

A8 

Estimated  Standard  Deviation  of  Forecast 
Total  Cost 

of 

!)7 

A9 

Values  for  (s^,S^) 

08  to  013 

A10 

Standard  Deviations  of  (s.,S^)  Values 

014  to  01 6 

Note:  For  corresponding  data  in  MacCormick  (1974),  see  tables 
of  the  same  number  in  Appendices  A of  those  reports. 
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Appendix  E 


Single-item  Data 


Model  IV 


Summary  of  Data  for  16  Items  with  Negative  Binomial  Demand 


Distributions  (Variance/Mean  = 3)  Controlled  with: 
Optimal  Policies  (DP) 

Power  Approximation  (PA) 

Statistical  Power  Approximation 


Table  A1  Average  Total  Cost 

A3  Period-End  Inventory 

A4  Period-End  Backlog  as  Proportion  of  Mean  Demand 


A5  Frequency  of  Periods  with  Backlog 
A6  Replenishment  Frequency 

A7  Estimated  Bias  of  Forecast  of  Total  Cost 


A8  Estimated  Standard  Deviation  of  Forecast  of 
Total  Cost 


A9  Values  for  (s^«Sj) 

A10  Standard  Deviations  of  (s^.S^)  Values 


E8  lo  El  3 


i:i4  to  i : 1 5 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  tables 
of  the  same  number  in  Appendices  A of  those  reports. 


Appendix  F 


Multi-item  Data  for  the  Power  Approximation 


Stationary  Model 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Optimally  with  Full  Information. 


Table  FI  Sources  of  Expected  Total  Cost  FI 

F3  Sources  of  Aggregate  Period-End  Inventory  F2 

F4  Sources  of  Aggregate  Backlog  Cost  F3 

F5  Sources  of  Aggregate  Replenishment  Cost  F4 

F6  Backlog  Frequency  F5 

F7  Weighted  Proportion  of  Demand  Backlogged  F6 

F8  Replenishment  Frequency  F7 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table  of  the 
same  number  in  his  Appendix  C. 


Appendix  G 


Multi-item  Data  for  the  Power  Approximation 
Stationary  Model 

System  of  16  Items  with  Negative  Binomial  Demand  Distribut ions 
(Variance/Mean  = 3)  Controlled  Approximately  Optimally  with 
Full  Information. 

page 


Table  G1  Sources  of  Expected  Total  Cost  G1 

G3  Sources  of  Aggregate  Period-End  Inventory  G2 

G4  Sources  of  Aggregate  Backlog  Cost  G3 

G5  Sources  of  Aggregate  Replenishment  Cost  G4 

G 6 Backlog  Frequency  G5 

G7  Weighted  Proportion  of  Demand  Backlogged  G6 

G8  Replenishment  Frequency  G7 

G15  Sources  of  Total  Cost  (%  Excess  Over  DP)  G8 

G16  Sources  of  Aggregate  Period-End  Inventory(%  Excess  Over  DP)  C9 

G17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  G10 

G18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Gil 

G19  Backlog  Frequency  (%  Excess  Over  DP)  C12 

G20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess  Over  DP)  G]3 

G21  Replenishment  Frequency  (%  Excess  Over  DP)  C14 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 

of  the  same  number  in  his  Appendices  E and  F. 


Appendix  H 

Multi-item  Data  for  the  (24,24)  Statistical  Power  Approximation 

Stationary  Model 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  with  Statistical  Information  from 
a 24-Period  Demand  History,  with  Revision  Every  24  Periods, 

Using  Regression  Estimates  of  Demand  Means  and  Variances. 


page 


Table  HI  Sources  of  Total  Cost  HI 

H3  Sources  of  Aggregate  Period-End  Inventory  H2 

H4  Sources  of  Aggregate  Backlog  Cost  H3 

H5  Sources  of  Aggregate  Replenishment  Cost  H4 

H6  Backlog  Frequency  H5 

H7  Weighted  Proportion  of  Demand  Backlogged  H6 

H8  Replenishment  Frequency  H7 

H15  Sources  of  Total  Cost  (%  Excess  Over  DP)  118 

H16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  H9 

H17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  H10 

H18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess 

Over  DP)  Hll 

H19  Backlog  Frequency  (%  Excess  Over  DP)  H12 

H20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  HI 3 

H21  Replenishment  Frequency  (%  Excess  Over  DP)  1114 


Note: 


For  corresponding  data  in  MacCormick 
of  the  same  number  in  his  Appendices 


(1974),  see  the  table 
K.  and  L. 


Appendix  I 


Stationary  Model 

Forecasting  Properties  of  Inventory  System  of  16  Items  with 
Negative  Binomial  Demand  Distributions  (Variance/Mean  = 3) 

Controlled  with  Statistical  Information,  Revision  Taking  Place 
Every  24  Periods  Using  a 24-Period  Demand  History  and  Regression 
Estimates  of  Demand  Means  and  Variances.  Forecasts  Made  at  Each 
Revision  Using  a 24-Period  Demand  History. 

£a£e 


Table  II  Sources  of  Forecast  of  Total  Cost  11 

12  Sources  of  Forecast  of  Aggregate  Period-End  Inventory  12 

13  Sources  of  Forecast  of  Aggregate  Backlog  Cost  13 

14  Sources  of  Forecast  of  Aggregate  Replenishment  Cost  14 

15  Forecast  of  Backlog  Frequency  15 

16  Forecast  of  Weighted  Proportion  of  Demand  Backlogged  16 

17  Forecast  of  Replenishment  Frequency  17 

18  Sources  of  Variance  of  Forecast  of  Total  Cost  18 

19  Sources  of  S.D.  of  Forecast  of  Total  Cost  19 

110  Sources  of  S.D.  of  Forecast  of  Aggregate 

Period-End  Inventory  110 

111  Sources  of  S.D.  of  Forecast  of  Aggregate  Backlog  Cost  111 

112  Sources  of  S.D.  of  Forecast  of  Aggregate  Replenishment  Cost  112 

113  S.D,  of  Forecast  of  Backlog  Frequency  113 

114  Sources  of  Total  Cost  (%  Underestimate  of  Actual  by 

Forecast)  1 14 

115  Sources  of  Aggregate  Period-End  Inventory  (%  Underestimate 

of  Actual  by  Forecast)  115 

116  Sources  of  Aggregate  Backlog  Cost  (%  Underestimate  of 

Actual  by  Forecast)  116 

117  Sources  of  Aggregate  Replenishment  Cost  (%  Underest  invite 

of  Actual  by  Forecast)  117 

118  Backlog  Frequency  (%  Underestimate  of  Actual  by  Forecast)  118 

119  Weighted  Proportion  of  Demand  Backlogged  (%  Underestimate 

of  Actual  by  Forecast)  119 

120  Replenishment  Frequency  (%  Underestimate  of  Actual  by 

Forecast)  120 


Note:  The  corresponding  appendices  In  MacCormick  (1974),  are  his 
Appendices  T and  U. 


Appendix  J 


Multi-item  Data  for  the  Power  Approximation 

Model  I 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Optimally  with  Full  Information. 


£i!ge 


Table  J1  Sources  of  Expected  Total  Cost  Jl 

J3  Sources  of  Aggregate  Period-End  Inventory  J2 

J4  Sources  of  Aggregate  Backlog  Cost  J3 

J5  Sources  of  Aggregate  Replenishment  Cost  J4 

J6  Backlog  Frequency  J5 

J7  Weighted  Proportion  of  Demand  Backlogged  J6 

J8  Replenishment  Frequency  J7 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendix  C. 
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Appendix  K 


Multi- Item  Data  for  the  Power  Approximation 


Model  I 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Approximately  Optimally  with  KulL 
Information. 

Table  K.1  Sources  of  Expected  Total  Cost  K.1 

K3  Sources  of  Aggregate  Period-End  Inventory  K2 

K4  Sources  of  Aggregate  Backlog  Cost  Kl 

K5  Sources  of  Aggregate  Replenishment  Cost  K4 

K6  Backlog  Frequency  K5 

K7  Weighted  Proportion  of  Demand  Backlogged  K6 

K8  Replenishment  Frequency  K7 

K15  Sources  of  Total  Cost  (%  Excess  Over  DP)  K8 

K16  Sources  of  Aggregate  Period-End  Inventory 

(%  Excess  Over  DP)  K9 

K17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  K10 

Kl8  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Kll 

K19  Backlog  Frequency  (%  Excess  Over  DP)  K12 

K20  Weighted  Proportion  of  Demand  Backlogged 

(%  Excess  Over  DP)  Ki  3 

K21  Replenishment  Frequency  (%  Excess  Over  DP)  K14 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  E and  F. 
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Appendix  L 

Multi-Item  Data  for  the  (24,24)  Statistical  Power  Approximation 

Model  I 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
•(Variance/Mean  = 3)  Controlled  with  Statistical  Information  from 
a 24-Period  Demand  History,  with  Revision  Every  24  Periods, 

Using  Regression  Estimates  of  Demand  Means  and  Variances. 

Pa-KV 


Table  LI  Sources  of  Total  Cost  LI 

L3  Sources  of  Aggregate  Period-End  Inventory  L2 

L4  Sources  of  Aggregate  Backlog  Cost  L3 

L5  Sources  of  Aggregate  Replenishmnet  Cost  L4 

L6  Backlog  Frequency  1.5 

L7  Weighted  Proportion  of  Demand  Backlogged  L6 

L8  Replenishment  Frequency  L7 

L15  Sources  of  Total  Cost  (%  Excess  Over  DP)  L8 

L16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  L9 

L17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  L10 

L18  Sources  o'  Aggregate  Replenishment  Cost  (%  Excess 

Over  DP)  Lll 

L19  Backlog  Frequency  (%  Excess  Over  DP)  L12 

L20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

(Over  DP)  1.13 

L21  Replenishment  Frequency  (%  Excess  Over  DP)  1.14 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  K and  L. 


Appendix  M 


Multi-item  Forecasts  for  the  (24,24)  Statistical  Power  Approximation 

Model  I 

Forecasting  Properties  of  Inventory  System  of  16  Items  with 
Negative  Binomial  Demand  Distributions  (Variance/Mean  = 3) 
Controlled  with  Statistical  Information,  Revision  Taking  Place 
Every  24  Periods  Using  a 24-Period  Demand  History  and  Regression 
Est invites  of  Demand  Means  and  Variances.  Forecasts  Made  at  Each 
Revision  Using  a 24-Period  Demand  History. 
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Table  Ml  Sources  of  Forecast  of  Total  Cost  Mi 

M2  Sources  of  Forecast  of  Aggregate  Period-End  Inventory  M2 

M3  Sources  of  Forecast  of  Aggregate  Backlog  Cost  M3 

M4  Sources  of  Forecast  of  Aggregate  Replenishment  Cost  M4 

M5  Forecast  of  Backlog  Frequency  M5 

M6  Forecast  of  Weighted  Proportion  of  Demand  BAcklogged  M6 

M7  Forecast  of  Replenishment  Frequency  M7 

M8  Sources  of  Variance  of  Forecast  of  Total  Cost  MB 

M9  Sources  of  S.D.  of  Forecast  of  Total  Cost  M9 

M10  Sources  of  S.D.  of  Forecast  of  Aggregate 

Period-End  Inventory  M10 

Mil  Sources  of  S.D.  of  Forecast  of  Aggregate  Backlog  Cost  Mil 

M12  Sources  of  S.D.  of  Forecast  of  Aggregate 

Replenishment  Cost  M12 

M13  S.D.  of  Forecast  of  Backlog  Frequency  M13 

M14  Sources  of  Total  Cost  (%  Underestimate  of  Actual  by 

Forecast)  M14 

M15  Sources  of  Aggregate  Period-End  Inventory  (X  Underestimate 

of  Actual  by  Forecast)  M15 

M16  Sources  of  Aggregate  Backlog  Cost  (%  Underestimate  of 

Actual  by  Forecast)  M16 

M17  Sources  of  Aggregate  Replenishment  Cost  (%  Underestimate 

of  Actual  by  Forecast)  M17 

M18  Backlog  Frequency  (%  Underestimate  of  Actual  by  Forecast)  M18 

M19  Weighted  Proportion  of  Demand  Backlogged  (%  Underest  invite 

of  Actual  by  Forecast)  M19 

M20  Replenishment  Frequency  (%  Underestimate  of  Actual  by 

Forecast)  M20 


Note:  The  corresponding  appendices  in  MacCormick  (1974),  are  his 
Appendices  T and  U. 


Appendix  N 


Multi-item  Data  for  the  Power  Approximation 


Model  II 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Optimally  with  Full  Information. 
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Table  N1  Sources  of  Expected  Total  Cost  N1 

N3  Sources  of  Aggregate  Period-End  Inventory  N2 

N4  Sources  of  Aggregate  Backlog  Cost  N3 

N5  Sources  of  Aggregate  Replenishment  Cost  N4 

N6  Backlog  Frequency  N5 

N7  Weighted  Proportion  of  Demand  Backlogged  N6 

N8  Replenishment  Frequency  N7 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendix  C. 


Appendix  0 


Multi-item  Data  for  the  Power  Approximation 


Model  II 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Approximately  Optimally  with 
Full  Information. 

Table  01  Sources  of  Expected  Total  Cost  01 

03  Sources  of  Aggregate  Period-End  Inventory  02 

04  Sources  of  Aggregate  Backlog  Cost  03 

05  Sources  of  Aggregate  Replenishment  Cost  04 

06  Backlog  Frequency  05 

07  Weighted  Proportion  of  Demand  Backlogged  06 

08  Replenishment  Frequency  07 

015  Sources  of  Total  Cost  (%  Excess  Over  DP)  08 

016  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  09 

017  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  010 

018  Sources  of  Aggregate  Replenishment  Cost  (%  Excess 

Over  DP)  Oil 

019  Backlog  Frequency  (%  Excess  Over  DP)  012 

020  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  013 

021  Replenishment  Frequency  (%  Excess  Over  DP)  014 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  E and  F. 


Appendix  P 


Model  II 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  with  Statistical  Information  from 
a 24-Period  Demand  History,  with  Revision  Every  24  Periods, 

Using  Regression  Estimates  of  Demand  Means  and  Variances. 
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Table  PI  Sources  of  Total  Cost  PI 

P3  Sources  of  Aggregate  Period-End  Inventory  P2 

P4  Sources  of  Aggregate  Backlog  Cost  P3 

P5  Sources  of  Aggregate  Replenishment  Cost  P4 

P6  Backlog  Frequency  P5 

P7  Weighted  Proportion  of  Demand  Backlogged  P6 

P8  Replenishment  Frequency  P7 

P15  Sources  of  Total  Cost  (%  Excess  Over  DP)  P8 

P16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  P9 

P17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  P10 

P18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Pll 

P19  Backlog  Frequency  (%  Excess  Over  DP)  P12 

P20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  PI 3 

P21  Replenishment  Frequency  (%  Excess  Over  DP)  P14 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 

of  the  same  number  in  his  Appendices  K and  L. 


Appendix  Q 


Multi-Item  Forecasts  for  the  (24, 24)  Statistical  Power  Approximation 

Model  II 

Forecasting  Properties  of  Inventory  System  of  16  Items  with 
Negative  Binomial  Demand  Distributions  (Variance/Mean  = 3) 

Controlled  with  Statistical  Information,  Revision  Taking  Place 
Every  24  Periods  Using  a 24-Period  Demand  History  and  Regression 
Estimates  of  Demand  Means  and  Variances.  Forecasts  Made  at 
Each  Revision  Using  a 24-Perlod  Demand  History. 
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Table  Q1  Sources  of  Forecast  of  Total  Cost  Q1 

Q2  Sources  of  Forecast  of  Aggregate  Period-End  Inventory  Q2 

Q3  Sources  of  Forecast  of  Aggregate  Backlog  Cost  Q3 

Q4  Sources  of  Forecast  of  Aggregate  Replenishment  Cost  Q4 

Q5  Forecast  of  Backlog  Frequency  Q5 

Q6  Forecast  of  Weighted  Proportion  of  Demand  Backlogged  Q6 

Q7  Forecast  of  Replenishment  Frequency  Q7 

Q8  Sources  of  Variance  of  Forecast  of  Total  Cost  Q8 

Q9  Sources  of  S.D.  of  Forecast  of  Total  Cost  Q9 

Q10  Sources  of  S.D.  of  Forecast  of  Aggregate 

Period-End  Inventory  Q10 

Qll  Sources  of  S.D.  of  Forecast  of  Aggregate  Backlog  Cost  Qll 

Q12  Sources  of  S.D.  of  Forecast  of  Aggregate 

Replenishment  Cost  <112 

Q13  S.D.  of  Forecast  of  Backlog  Frequency  Q13 

Q14  Sources  of  Total  Cost  (%  Underestimate  of  Actual  by 

Forecast)  Q14 

Q15  Sources  of  Aggregate  Period-End  Inventory  (%  Underestimate 

of  Actual  by  Forecast)  Q15 

Q16  Sources  of  Aggregate  Backlog  Cost  (%  Underestimate  of 

Actual  by  Forecast)  i)lb 

Q17  Sources  of  Aggregate  Replenishment  Cost  (%  Underest  invite 

of  Actual  by  Forecast)  Q17 

Q18  Backlog  Frequency  (%  Underestimate  of  Actual  by  Forecast)  Q18 

Q19  Weighted  Proportion  of  Demand  Backlogged  (%  Underestimate 

of  Actual  by  Forecast)  Q19 

Q20  Replenishment  Frequency  (%  Underestimate  of  Actual  by 

Forecast)  Q20 


Note:  The  corresponding  appendices  in  MaeCormiek  (1974),  are  his 

Appendices  T and  U. 


Appendix  R 


Multi-item  Data  for  the  Power  Approximation 
Model  III 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Optimally  with  Full  Information. 


pa£e 


Table  R1  Sources  of  Expected  Total  Cost  R1 

R3  Sources  of  Aggregate  Period-End  Inventory  R2 

R4  Sources  of  Aggregate  Backlog  Cost  R3 

R5  Sources  of  Aggregate  Replenishment  Cost  R4 

R6  Backlog  Frequency  R5 

R7  Weighted  Proportion  of  Demand  Backlogged  R6 

R8  Replenishment  Frequency  R7 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendix  C. 
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Appendix  S 


Multi-item  Data  for  the  Power  Approximation 

Model  HI 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Approximately  Optimally  with 
Full  Information. 

page 


Table  SI  Sources  of  Expected  Total  Cost  SI 

53  Sources  of  Aggregate  Period-End  Inventory  S2 

54  Sources  of  Aggregate  Backlog  Cost  S3 

S3  Sources  of  Aggregate  Replenishment  Cost  S4 

56  Backlog  Frequency  S5 

57  Weighted  Proportion  of  Demand  Backlogged  S6 

58  Replenishment  Frequency  S7 

S13  Sources  of  Total  Cost  (%  Excess  Over  DP)  SR 

516  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  UP)  S9 

517  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  S10 

SIR  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Sll 

519  Backlog  Frequency  ( X Excess  Over  DP)  S12 

520  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  S13 

521  Replenishment  Frequency  (%  Excess  Over  DP)  SI 4 


Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  E and  F. 


Appendix  T 


Multi-item  Data  for  the  (24,24)  Statistical  Power  Approximation 

Model  III 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  with  Statistical  Information  from 
a 24-Period  Demand  History,  with  Revision  Every  24  Periods, 

Using  Regression  Estimates  of  Demand  Means  and  Variances. 

Table  T1  Sources  of  Total  Cost  T1 

T3  Sources  of  Aggregate  Period-End  Inventory  T2 

T4  Sources  of  Aggregate  Backlog  Cost  T3 

T5  Sources  of  Aggregate  Replenishment  Cost  T4 

T6  Backlog  Frequency  T5 

T7  Weighted  Proportion  of  Demand  Backlogged  T6 

T8  Replenishment  Frequency  T7 

T15  Sources  of  Total  Cost  (%  Excess  Over  DP)  T8 

T16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  T9 

T17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  TIO 

T18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Til 

T19  Backlog  Frequency  (%  Excess  Over  DP)  T12 

T20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  T13 

T21  Replenishment  Frequency  (%  Excess  Over  DP)  T14 


Note:  For  corresponding  data  in  MacCormlck  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  K and  L. 


Appendix  U 


Multi-item  Forecasts  for  the  (24,24)  Statistical  Power  Approxlma t ion 

Model  III 

Forecasting  Properties  of  Inventory  System  of  16  Items  wiLh 
Negative  Binomial  Demand  Distributions  (Variance/Mean  = 3) 
Controlled  with  Statistical  Information,  Revision  Taking  Place 
Every  24  Periods  Using  a 24-Period  Demand  History  and  Regres- 
sion Estimates  of  Demand  Means  and  Variances.  Forecasts  Made 
at  Each  Revision  Using  a 24-Period  Demand  History. 
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Table  U1  Sources  of  Forecast  of  Total  Cost  U1 

U2  Sources  of  Forecast  of  Aggregate  Period-End  Inventory  U2 

U3  Sources  of  Forecast  of  Aggregate  Backlog  Cost  U3 

U4  Sources  of  Forecast  of  Aggregate  Replenishment  Cost  IJ4 

U5  Forecast  of  Backlog  Frequency  U5 

U6  Forecast  of  Weighted  Proportion  of  Demand  Backlogged  U6 

U7  Forecast  of  Replenishment  Frequency  U7 

U8  Sources  of  Variance  of  Forecast  of  Total  Cost  U8 

U9  Sources  of  S.D.  of  Forecast  of  Total  Cost  U9 

U10  Sources  of  S.D.  of  Forecast  of  Aggregate  Period-End 

Inventory  U10 

Ull  Sources  of  S.D.  of  Forecast  of  Aggregate  Backlog  Cost  Ull 

U12  Sources  of  S.D.  of  Forecast  of  Aggregate  Replenishment  Cost  U12 

U13  S.D.  of  Forecast  of  Backlog  Frequency  U13 

U14  Sources  of  Total  Cost  (%  Underestimate  of  Actual  by 

Forecast)  U14 

U15  Sources  of  Aggregate  Period-End  Inventory  (X  Underestimate 

of  Actual  by  Forecast)  U15 

U16  Sources  of  Aggregate  Backlog  Cost  (X  Underestimate  of 

Actual  by  Forecast)  U16 

U17  Sources  of  Aggregate  Replenishment  Cost  (X  Underest  invite 

of  Actual  by  Forecast)  U17 

L'18  Backlog  Frequency  (X  Underestimate  of  Actual  by  Forecast)  U18 

U19  Weighted  Proportion  of  Demand  Backlogged  (X  Underestimate 

of  Actual  by  Forecast)  U19 

U20  Rep lenishment  Frequency  (X  Underestimate  of  ActuaL  by 

Forecast)  U20 


Note:  The  corresponding  appendices  in  MacCormick  (1974),  are  his 
Appendices  T and  U. 


Appendix  V 

Multi-item  Data  for  the  Power  Approximation 
Model  IV 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Optimally  with  Full  Information. 
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Table  VI 

Sources  of  Expected  Total  Cost 

VI 

V3 

Sources  if  Aggregate  Period-End  Inventory 

V2 

V4 

Sources  of  Aggregate  Backlog  Cost 

V3 

V5 

Sources  of  Aggregate  Replenishment  Cost 

V4 

V 6 

Backlog  Frequency 

V5 

V7 

Weighted  Proportion  of  Demand  Backlogged 

V6 

V8 

Replenishment  Frequency 

V7 

Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendix  C. 


Appendix  W 


Multi-item  Data  for  the  Power  Approximation 


Model  IV 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  Approximately  Optimally  with 
Full  Information. 
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Table  W1  Sources  of  F.xpected  Total  Cost  W1 

W3  Sources  of  Aggregate  Period-End  Inventory  W2 

W4  Sources  of  Aggregate  Backlog  Cost  W3 

W5  Sources  of  Aggregate  Replenishment  Cost  W4 

W6  Backlog  Frequency  W5 

W7  Weighted  Proportion  of  Demand  Backlogged  W6 

W8  Replenishment  Frequency  W7 

W15  Sources  of  Total  Cost  (%  Excess  Over  DP)  W8 

W16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess  W9 

Over  DP) 

W17  Sources  of  Aggregate  Backlog  Cost  ( [X  Excess  Over  DP)  W10 

W18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess  Over  DP)  Wll 
W19  Backlog  Frequency  (%  Excess  Over  DP)  W12 

W20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess  Over  DP)  W13 

W21  Replenishment  Frequency  (%  Excess  Over  DP)  W14 


Note:  For  correspond! ng  data  in  MacCormick  (1974),  see  the  table  of 

the  same  number  in  his  Appendices  E and  F. 


Appendix  X 


Multi-item  Data  for  the  (24,24)  Statistical  Power  Approximate  ion 

Model  IV 

System  of  16  Items  with  Negative  Binomial  Demand  Distributions 
(Variance/Mean  = 3)  Controlled  with  Statistical  Information  from 
a 24-Period  Demand  History,  with  Revision  Every  24  Periods, 

Using  Regression  Estimates  of  Demand  Means  and  Variances. 
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Table  XI  Sources  of  Total  Cost  XI 

X3  Sources  of  Aggregate  Period-End  Inventory  X2 

X4  Sources  of  Aggregate  Backlog  Cost  X3 

X5  Sources  of  Aggregate  Replenishment  Cost  X4 

X6  Backlog  Frequency  X5 

X7  Weighted  Proportion  of  Demand  Backlogged  X6 

X8  Replenishment  Frequency  X7 

X15  Sources  of  Total  Cost  (%  Excess  Over  DP)  X8 

X16  Sources  of  Aggregate  Period-End  Inventory  (%  Excess 

Over  DP)  X9 

X17  Sources  of  Aggregate  Backlog  Cost  (%  Excess  Over  DP)  X10 

X18  Sources  of  Aggregate  Replenishment  Cost  (%  Excess 

Over  DP)  Xll 

X19  Backlog  Frequency  (%  Excess  Over  DP)  X12 

X20  Weighted  Proportion  of  Demand  Backlogged  (%  Excess 

Over  DP)  X13 

X21  Replenishment  Frequency  (%  Excess  Over  DP)  X14 

Note:  For  corresponding  data  in  MacCormick  (1974),  see  the  table 
of  the  same  number  in  his  Appendices  K and  L. 


Appendix  Y 


Multi-item  Forecasts  for  the  (24,24)  Statistical  Power  Approximation 

Mod  el  IV 

Forecasting  Properties  of  Inventory  System  of  16  Items  with 
Negative  Binomial  Demand  Distributions  (Variance/Mean  = 3) 
Controlled  with  Statistical  Information,  Revision  Taking  Place 
Every  24  Periods  Using  a 24-Period  Demand  History  and  Regres- 
sion Estimates  of  Demand  Means  and  Variances.  Forecasts  Made 
at  Each  Revision  Using  a 24-Period  Demand  History. 
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Table  Y1  Sources  of  Forecast  of  Total  Cost  Y1 

Y2  Sources  of  Forecast  of  Aggregate  Period-End  Inventory  Y2 

Y3  Sources  of  Forecast  of  Aggregate  Backlog  Cost  Y3 

Y4  Sources  of  Forecast  of  Aggregate  Replenishment  Cost  Y4 

Y5  Forecast  of  Backlog  Frequency  Y5 

Y6  Forecast  of  Weighted  Proportion  of  Demand  Backlogged  Y6 

Y7  Forecast  of  Replenishment  Frequency  Y7 

Y8  Sources  of  Variance  of  Forecast  of  Total  Cost  Y8 

Y9  Sources  of  S.D.  of  Forecast  of  Total  Cost  Y9 

Y10  Sources  of  S.D.  of  Forecast  of  Aggregate  Period-End 

Inventory  Y10 

Yll  Sources  of  S.D.  of  Forecast  of  Aggregate  Backlog  Cost  Yll 

Y12  Sources  of  S.D.  of  Forecast  of  Aggregate 

Replenishment  Cost  Y12 

Y13  S.D.  of  Forecast  of  Backlog  Frequency  Y13 

Y14  Sources  of  Total  Cost  (%  Underestimate  of  Actual  by 

Forecast)  Y14 

Y15  Sources  of  Aggregate  Period-End  Inventory  (%  Underestimate 

of  Actual  by  Forecast)  Y15 

Y16  Sources  of  Aggregate  Backlog  Cost  (%  Underestimate  of 

Actual  by  Forecast)  Y16 

Y17  Sources  of  Aggregate  Replenishment  Cost  (%  Underestimate 

of  Actual  by  Forecast)  Y17 

Y18  Backlog  Frequency  (%  Underestimate  of  Actual  by  Forecast)  Y18 

Y19  Weighted  Proportion  of  Demand  Backlogged  (%  Underestimate 

of  Actual  by  Forecast)  Y19 

Y20  Replenishment  Frequency  (%  Underestimate  of  Actual  by 

Forecast)  Y20 


Note:  The  corresponding  appendices  in  MacCormick  (1974),  are  his 
appendices  T and  U. 


